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Doctor of Science. : 


ABoOTRACT 


The problem studied in this thesis is the guidance of interplane- 
tary vehicles which are thrusting for a large portion of the transfer. 
The vehicle is represented by a seven component state vector Consiga 
ing of the position, velocity and mass of the spacecraft. The analysis 
is linearized by assuming that the actual state of the vehicle differs 
only a small amount from a known reference state. The reference 
trajectory is assumed to be a propellant-optimal path connecting the 
initial and final points. 


The goal of the postulated guidance system is to satisfy position 
and velocity conditions at the target with minimum propellant expedi- 
ture. Both fixed-time-of-arrival and variable-time-of-arrival guides 
ance are discussed. Specification of the guidance criterion taiime 
above manner permits the techniques of optimal control theory to be 
applied to the problem. Emphasis is placed on finding an analytic 
solution of the linearized equations. The desired solution is the con- 
trol program which satisfies boundary conditions and minimizes pro- 
pellant expenditure. 


The method for solving the guidance problem is shown to be suit- 
able as a technique for computing optimal reference trajectories. The 
trajectories are computed by iterative application of the guidance sol- 
ution. Application of the guidance solution to the trajectory problem 
is shown to exploit an interpretation of the Euler equations which per- 
mits simplification of the computation technique. 


The guidance solution is tested in a numerical example by using 
it to compute trajectories from Earth's orbit to the Martian orbit for 
different low-thrust vehicles. 


The guidance solution is based on the assumption that vehicle 
state is known at the time a new control program is to be generated. 
Prior studies by several investigators detail methods of using celestial 
measurements to estimate state. A portion of this report is devoted 
to extending the method of celestial measurements to include measure- 
ment of engine performance. The additional measurement is shown 
to improve the estimate of state. 
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A discussion is presented of the difficulties arising from differ- 
milecs il tae Crirerion lor Optimality as interpreted from the calculus 
of variations and from Pontryagin's maximum principle. 
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CHAPTER I 


INTRODUCTION 


Early in the investigation of rocket propulsion for extra-terrestrial 
imeavel, it became apparent that conventional chemical propellanis were 
inadequate for many interesting space missions Bere hits fact is 
due to the relatively low energy content per unit mass (specific energy) 
Sieewemical fuels, The theoretical mass ratios required for many inter- 
eels missions in the solar system approach numbers of the order of 


10° and higher when chemical fuels are used . 


ime need tor more efficient means of propulsion has led to the 
investigation of energy sources other than chemical reactions. Such 
Ewuenes have produced an entire spectrum of propulsive techniques®, 
each having particular mission capabilities and each having its own 
Mmeoretical and practical difficulties. In general, the use of higher 
Beecilic energies is accompanied by a decrease in propellant flow rate 
and longer propulsion time for a given mass ratio. In limiting cases 
the propulsion time equals the transfer time. In addition, a longer pro- 
Palesian time permits reduction of thrust levels for a given total impulse. 
The terms "continuous thrust’ and "low thrust'' stem from these two 
elfects of high specific energy. Although the terms are not synonymous, 
they are often used interchangeably in the literature since they apply to 
the same types of vehicles. In the subsequent discussion, “continuous 
thrust'' is used to describe vehicles which thrust for a major portion of 
the transfer and ‘low thrust'' refers to vehicles which have acceleration 


levels less than about 10°°g. 


This study is concerned primarily with those propulsion methods 
which rely upon a separate energy source* for the generation of pro- 


pulsion energy. With few exceptions, these are low thrust devices. 


“Underlined words are defined in footnotes. 


SCPalele ener sy source = preopulsiOl cy siem ellaractetizea oy 2 power 
plant which is independent of the thrust-producing mechanism. 


Where it is necessary to be more specific, anon-board nuclear reactor 


is assumed. 


We shall not be concerned with further justifying the use of these 
devices, nor with defining their regions of usefulness. Such questions 


are well covered in the intevatunes 


The goal here is to assume 
that such vehicles will exist and to study the problem of guidance ir- 


respective of the mission or of the particular propulsion method used. 
iol) Ghne Guidance. .oplen 


The presence of a thrust acceleration, acting over a significant 
portion of the spacecrait trajectory, introduces complexities in sea 
analysis of spacecraft motion which make many of the techniques in 
common use for ballistic vehicles inapplicable. In particular, theveves 
gant conic representations can be used only when the thrust may be 
treated as a small disturbing force. In general this is the case only in 
the region near a central body. Consequently it is desirable to fornia 
late continuous-thrust investigations within a body of mathematics hav- 
ing sufficient applicability to treat most problems of interest, Tae 
calculus of variations and the concepts of optimal control theory meet 


this need. 


In the terminology of this branch of applied mathematics, guidance 
of continuous-thrust spacecraft is the problem of finding a control pro- 
gram which transfers the spacecraft between two given points in space 
subject to constraints on the velocity at both terminal points. Since 
such programs are not unique, it is possible to place additional con- 
straints on the solution such that some desired quantity (the cost func- 
tion) is maximized or minimized. Thus, in addition to satisfying fixed 
boundary conditions at both launch and target points, solutions may be 


found which minimize transfer time, or propellant consumed or some 


guidance: as used in this thesis, the term refers to the process of 
determining a control program. 


control: as used in this thesis, the term refers to the quantity or quan- 
tities (thrust or acceleration) representing the propulsive effort. It 
also refers to the mechanical process of directing the propulsive effort 
in accordance with the control program. 


other quantity appropriate to the mission This thesis is concerned 
with maximizing only the final mass. At the present time no method 
has been found which will prove rigorously that the solutions obtained 
in inverse square (and more general) gravitational fields satisfy the 
condition for an absolute maximum of the final mass’, One can only 
show that the solutions produce a maximum in the region of space under 
consideration, i.e. local maxima, and then demonstrate from physical 


reasoning that other solutions are not likely. 


From the viewpoint of guidance, initial conditions for the two-point 
boundary value problem are represented by the present state of the 
vehicle. Ifa vehicle can be flown in accordance with the solution that 
exists at the instant of launch, only one solution is needed Such a case 
is a problem in control, not in guidance. That is, the desired thrust 
program is known (guidance) and the magnitude and direction of thrust 
must be controlled so that the spacecraft follows the desired trajectory 
(control). Even with ‘tight’ control locps. however, perturbations in- 
evitably occur which cause the spacecraft to depart from tne original 
optimized trajectory. When this happens a new control program must 
be found which will cause transfer from the present state to the final 
eae such that final mass 1s Maximized. This is the problem we desire 
to solve. it follows, that if a solution can be found for the guidance prob- 
lem, then by considering the launch state as the present state, the op- 
timized trajectory connecting the launch point and target point may be 
determined This latter application of the guidance techniques leads to 


that part of the thesis described as ‘trajectory computation " 


In the preceding paragraph it was emphasized that the present state 
of the vehicle serves as the initial condition for the boundary value prob- 
lem It is necessary, therefore, that some method exist for determining 


the state of the vehicle at any time. In this thesis the unique features 


slate, relers to the seven quantities which describe the vehicle in terms 
of physical variables 


present state: reiers to the vehicle state as determined by the space 
navigator, 1 e. the state at the present instant of time In this thesis 
navigation refers to the process of determining state from measure- 
ments. 


of a continuous-thrust vehicle are examined from the viewpoint of 


estimating vehicle state. 


Finally, the concept of constant exhaust power is examined as a 
criterion for low-thrust transfers. It is easily shown that if a linear 
system is power-limited, as in the case of separately powered rockets, 
the minimum expenditure of energy results from operation at maximum 
continuous power. The application of this principle has been used ex- 


679,10 


tensively tarousnouline diierauiiT. An analogous treatment of 


thrust-limited rockets 1s examined here. 
le2 Priore Studies 


Published works on guidance of interplanetary low-thrust vehicles 
were alimost NONECXIStea) pri@ra oscar lvelooo. Teller producedya 
guidance technique for cis-lunar space in his doctoral research in 1961. 
This technique consists of spiralling out from the Earth to some point 
from which the vehicle may coast to the vicinity of the moon, then 
matching the unique velocity vector, corresponding to the present poom. 
tion of the vehicle, which will result in achieving the target point. 
Miller showed that guidance of this type results ina relatively small 


fuel penalty for lunar missions. 


Friedlander’! formulated the problem in the classical calculus of 
variations and solved the adjoint equations in two dimensions for the 
sensitivity coefficients of the state variables along an optimized tra- 
jectory to Mars. In reference 12, he sulesesis a linearized coliiien 
which minimizes a quadratic function of the control variable variation. 
The solution approaches but does not attain the control program for 
maximum final mass. In reference 13 Friedlander applies his tech- 


niques to a vehicle using a Snap-8 power source. 


Ihe most recent work is thar of Pfeitier who applies some of the 
newer developments in the theory of optimal control to the low-thrust 
guidance problem. Pfeiffer solves the guidance problem and minimizes 
power limited: refers to a propulsion system characterized primarily 
by a maximum power level. 


thrust limited: refers to a propulsion system characterized primarily 
by a maximum thrust level 


a penalty function which is ‘equivalent’ to a quadratic form of the final 
state error. The control program produced by this method satisfies 
the boundary conditions ''as closely as possible using the penalty func- 


'' (the quotation marks are from the reference) Pfeiffer's method, 


tion. 
however, is not applicable to problems where certain boundary values 


are fixed. 


MlechinOr the recent work in Optimal control theory concerns sys- 
tems which have mathematical models similar to those for low-thrust 
vehicles. Many of the ideas developed in these investigations are 
Sumeciiy applicable to the problem considered here. Significant con- 
tributions are attributable to Pontryagin!?, Kalman?°® and Breakwell’ 
who have formulated existence theorems and derived necessary and 
sufficient conditions for optimal trajectories. Breakwell has also 
contributed important work in specifying the form of solutions with 
constrained control vectors Athans, Falb and fhacGrosa:” working 
together and individually have solved many special cases of optimal 


Beejecrories for constrained control vectors. 


Fundamental to any guidance study is a qualitative knowledge of the 
trajectory along which the space vehicle is to travel. Quite properly 
then, the earliest work in low-thrust propulsion consisted of studies of 
emaeide Characteristics and trajectory characteristics Several of the 
famlier studies of engine characteristics have been previously refer- 
eaced. IO those must be added the contributions of Langmuir’? and 


20 
Irving 


In the area of trajectory studies Tsien*! Detlor Med seme On sunc 
earliest work (1953). This was followed several years later with con- 
tributions by Lawden““, Moeckel?”, Melbourne“ “and Zimmerman, 
McKay and Rossa” The problem which confronted these authors was 
that analytic solutions to the trajectory problem can be found only for 
linear gravitational fields. For central force fields and more complex 
configurations, solving the two-point boundary value problem was a 
teqious trial-and-error procedure requiring the use of high speed digit - 
al computation. Satisfying an additional constraint for an optimized 


trajectory was tedious and time consuming even with high speed 


Computers, | Lhe worker Bryson?° and others in the late 1350"s and 
early 1960's served to simplify the machine procedures so that com- 
putation of optimal trajectories became less tedious and less time 


consuming. 
1.3 Thesis Philosophy and the Method of Approach 


Throughout the research and writing of this thesis the author has 
attempted to consider the low-thrust guidance problem from the view- 
point of a space navigator who is responsible for the safe and timely 
arrival of the spacecrait at the target point. The exirapolation of air= 
craft navigation experience into space navigation is, at best, a hazard- 
ous undertaking; however, it does provide a basis for certain decisions 
which have influenced the author's approach to this investigation. The 
following criteria were established from this philosophy and have been 


used when it became necessary to make definite assumptions. 


1) The mission is manned, probably utilizing more than one 


vehicle, each of which is manned by several crewmen. 


2) The mission duration is limited by consideration of human 


tolerances. 


3) The spacecraft configuration and the mission have been speci- 
fied. Hopefully, the spacecrait characteristics are Optimum 


for the mission, but may not be. 


4) Whatever the mission, at each of several points along the tra- 

jJeclory the Navigator dias ir cercholce.: 

a) to rendezvous with the target point at the preplanned time 
such that the trajectory minimizes propellant consumed, 

b) to rendezvous with the target point utilizing a time and 
trajectory such that propellant consumption is minimized. 

c) to rendezvous with the target point in minimum time Ute 
ing the available propellant. This alternative is not 


treated in the thesis. 


These guidelines establish the general context within which the guidance 
problem is to be solved. Let us now proceed to examine the specific 


items which complicate the solution. 


The mathematical theory concerning optimal trajectories dic- 
res tiat the iirst variation of the optimized quantily must 
vanish for the optimal path when the control is unconstrain- 
eg? (228,29 


state variable, the matrix of coefficients relating the control 


As a consequence, if the optimized quantity is a 


variables to the state variables is singular and its time inte- 
12,14 The 


Singularity is proven in Appendix C. Solutions for the optimal 


gral along the optimal trajectory is also singular 


conmeolwcvally require inversion Of This matrix. Phe prep 
lem of singular matrices is handled in this study by a method 
of deleting certain matrix elements which create the singular - 
ity and by the formation of a new matrix which can be inverted. 
The deletion method 1s an important part of the thesis and 
provides a general method for treating certain singularities 


without reformulating the problem. 


Optimal trajectories for constrained control variables often 

possess discontinuities in the first derivatives of one or more 
State variables and in the control variables. This mathemati- 
cal problem is handled by the use of switching functions which 


are continuous. 


Optimal trajectories for constrained control variables require 
Petlocs Of iInaximum Control magniiude, “Theretore, it ie 
maximum propulsive effort 1s required for the optimum ira- 
jectory, guidance around the optimum is limited to changes in 
thrust direction unless reserve propulsive power is available 
for guidance. The assumption that reserve power is available 


Peevceanorwihis sic” 


It was stated in section 1 1 that the low-thrust guidance problem 


Ol evariations 


ters. 


Conirol 


may be treated as a two-point boundary value problem in the calculus 


- (528,29 Thus techniques of the calculus of variations 


constitute a primary mathematical tool One of these techniques is the 
method of adjoints, which plays a fundamental part in subsequent chap- 


The author has borrowed heavily from newer theories in optimal 
Ligh 


Since the state space formulations widely used in the 


literature of that field are applicable to low-thrist euidances Ones, 
the more useful tools in optimal control theory is associated with 


Pontryagin’ a although other authors have used the same principle™ ae 


A derivation of Pontryagin's maximum principle is outlined in 


Appendix D for convenience of the reader. 


To facilitate notation and to preclude the possibility of fundamental 
notions becoming obscured by the quantity of algebraic detail, matrix 
notation and the ideas of matrix calculus are used throughout. Finally, 


to test the thesis, numerical analysis and an IBM 7094 were employed. 
1.4 Relationship to Prior Studies 


In the research preceding this study, the author began with the 
formulations of Friedlander!® and Melbourne’, and attempted to 
extend their ideas into areas of more general application and to find 
solutions to the guidance problem which were useful from the space 
navigator's viewpoint. One of the fundamental considerations was to 
find control laws which optimize propellant consumption. Cost func- 


tions which produce near-optimal propellant consumption were rejected. 


The idealized formulations for the separately powered rocket re- 
quire a wide range of thrust and of specific impulse as the vehicle 
traverses its trajectory. Curreni technology indicates that variable- 
specific-impulse thrusters will not be available in the forseeable 
future, at least for electrostatic vehicles. To Satisfy this engincermm 
restriction, investigators have continued to use the power-limited for- 
mulations but approximate the optimal thrust magnitude programs with 
regions of constant specific impulse where relatively low values of 
specific impulse are required and with coast elsewhere. It is shown 
in this thesis that the so-called 'bang-bang" control used to satisfy the 
engineering restriction can be derived by abandoning the concept of 
power-limited thrusting. 


In section 1.1 the necessity of estimating the state is discussed. 


Te BVODKS' OF een oes and Fetter and eters serve as the 


electrostatic (propulsion): A propulsion method depending upon the 
acceleration of charged particles through an electrostatic field. 


elertine point for exiending, to the low-thrust case, the navigation tech- 


niques (estimate of state) which have been developed for ballistic vehicles. 


Finally, using the guidance solutions an iterative technique for com- 
puting low-thrust trajectories is developed. Recent investigations in the 
field of trajectory computation appear to rely heavily on steepest ascent 
feenmmiques. These techniques require rather complex programs, care- 
imenanaling and suifer from the problem of slow convergence as the 
Serial trajectory is approached. The technique suggested in this 


study appears to be faster even for discontinuous control variables. 
im@eoummary of Thesis Objectives 


In closing this introductory chapter it seems appropriate to provide 
a sketch of the objectives of subsequent discussions. Briefly, the ob- 
Weeulves Oi this thesis are: 1) to present a linear. noniterative method 
ieeeaeriving a control program for guidance of low-thrust space vehicles 
with respect to a propellant-optimal reference, 2) to show that the com- 
puted control programs satisfy the necessary conditions for an optimum 
3) to show that the method may be extended to trajectory computation by 
successive iteration of the guidance solution; 4) to examine a method of 
estimating the state of continuous: thrusi vehicles and 5) to examine the 
concepts of power-limited and thrust-limited vehicles from the view- 


point of guidance. 


PicmuMmdamental arsyuiment Of the thesis nay be exiracted from 
these several objectives. It is: ‘'There exists a linear method which 
produces a propellant-optimal control program in a noniterative form 
for guidance of power-limited and thrust-limited space vehicles, and 
which provides a simple, rapidly converging iterative technique for com- 


puting propellant-optimal trajectories. | 


propellant-optimal: refers to a trajectory which results in minimum 
propellant usage. 


CHAPTER ii 


THE LOW-THRUST SPACE VEC rE 


2 t,oummary ol Chapter Il 


In this chapter the parameters which characterize low-thrust vehicle 


performance are derived and discussed. The constant-exhaust-power 
concept, which has been widely used in the past, is shown to be the op- 
timum method of engine control. The propellant cost accruing from 
engineering restrictions on variable specific impulse is computed tam 


field-free space by analytic methods. Both methods of engine control, 


the ideal and the practical, are discussed from the viewpoint of guidance. 


Finally, the results of a numerical study, which confirm the derivations, 


Bre presented, 
2.2 ne Use Oia Separate Eneroy-ouvec 


One may verify from momentum considerations that the instantane- 


ous force exerted on a space vehicle by its exhaust stream is: 
ie mc (25 


where f is the force vector, c is the oppositely directed exhaust velocity 


and m is the mass rate of the vehicle. (-m is the propellant flow rate. ) 


pince, tne force magnitude may bel nel cece nauanue m and C else 
varied inversely, the selection of a high exhaust velocity and low pro- 
pellant flow tends to reduce the total amount of propellant required for 
a given impulse. Unfortunately, processes which use the products of 
combustion as the propellant are unable to produce, simultaneously, the 
low flow rates and high exhaust velocities required for many interesting 
missions. If we use specific impulse, Isp? aS a measure of the engine 


performance, where 


Se ere (2-2) 
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and go issttemacceleramonm due tO gravity at the Earth's surface, then 
chemical systems are limited to values of specific impulse under 600 

seconds and direct nuclear systems to values of about 1600 seconds” | 
However, if a separate energy source is available, the source energy 

miay be Converted into electrical energy and used to accelerate charged 

propellant particles. Values of specific impulse in the range of 3000 


to 20,000 seconds appear attainable in this way” 


several energy sources and conversion processes are under inves- 
tigation. These may be divided into the two broad categories of direct 
ieimcdirect energy conversion. Direct methods may be characterized 
by the absence of a mechanical phase in the conversion process. Power 
magn SOlar cells is one example of this type. Direct production of e!ec- 
trical energy in a nuclear reactor is another. A proposal for this latter 
fethod is discussed in reference 33 however it has not been proven in 


the laboratory. 


Picirect Conversion of the Scurce energy into electrical energy 
appears, at the present time, to be more realistic for large manned 
spacecraft which have power requirements in the range of several 
megawatts The process generally considered most promising uses a 
nuclear source to power turbomachinery for the production of electrical 


energy. A block diagram of this type system is presented in Figure 2-a 


In this report, an on-board energy source is explicitly assumed 
although much of the guidance analysis is independent of such considera- 
tions. This assumption permits the power availability to be independent 
of the trajectory. This is not the case, for example, when solar energy 
Bemleciors are used since for a given collector area the power availa- 


bility varies inversely with the square of the distance to the sun. 
2.3 The Constant-Power Concept 


separate energy sources are most often described as power -limited 
devices. That is, their rate of energy conversion is the design criterion. 
If the energy produced by the source is all converted to propulsive en- 


ergy then the power in the exhaust stream is given by 


ioe 
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Fig. 2-a. Schematic diagram of separately powered rocket. 


-mc* 
Poa oe (2-3) 
2 
(The efficiency of the conversion process will be introduced in a sub- 


sequent discussion. ) 


It is well known that for a large class of problems in linear fields, 
minimum total energy expenditure results from operation at maximum 
continuous power. This result may be obtained for conservative fields 
in general by use of variational techniques as shown in reference 7. The 
usual assumption in low-thrust investigations, that exhaust power is a 
constant and equal to its maximum value, is therefore, quite valid in the 
idealized problem. This result is used to establish parameters for the 


separately powered rocket. 


Combining equation (2-3) with the relations 


ie (2-4) 
f=ma (23) 
one obtains 
“aa 2 
Se ae (eG) 
2 
m 2p 


megere 2a is the thrust acceleration of the vehicle. Integration of equa- 
tion (2-6) over the flight path yields 


1 


MN 


ee «an at (e= 
f mM, Zp. <0 


where the subscripts indicate initial and final times. Equation (2-7) is 


more conveniently expressed in terms of the mass ratio. 
i 
m i 
noeee sa f a’ dt (2G) 
2p od 


Three parameters will now be defined 


te 
J = f arat pac 
0 
mY) 
a = 2 ee (2-10) 
p 
Fry 
Soe (2-11) 
mM 
O 


J 1s the well known acceleration integral, a is the specific mass of the 
Peepulsion system and must include the total efficiency of energy con- 
version when p is the desired exhaust power; B is the propulsive system 
mass fraction and my is the propulsive system mass. Substituting 
these parameters into equation (2-8) one obtains 


Ni? =e ee (ih 
B 
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ACCELERATION MAGNITUDE (x04g,) 


For a spacecraft with fixed parameters and a specified mission, 
the minimum mass ratio and thus the minimum propellant usage is 


achieved for trajectories which make J a minimum. 


The objective of trajectory computation and of guidance is to fm 


a path which accomplishes the mission and minimizes J, 
2.4 Limitation of the Constant-Power Concept 


Using the equations of section 2. 3, the thrust acceleration may be 
expressed as 


eee (2-13) 


Civ) 


Figure 2-b shows a plot of the thrust acceleration magnitude during a 


fast transfer to Mars for which J is a minimum. The absolute valige 
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Fig. 2-b. Optimum acceleration level for variable-specific-impulse Earth-Mars transfer. 
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of acceleration has a variation over two orders of magnitude. From 
equation (2-13) it is apparent that the exhaust velocity, (i.e., the specif- 
ic impulse) must have a similar range of values since m is monotoni- 
cally decreasing. However, the current thinking of investigators who 
are studying propulsive devices indicates that in a given thrust pro- 
ducer, only a very small range of exhaust velocity variation will be 
feasible unless major advances in technology are forthcoming. The 
thrust producer can be designed for one particular exhaust velocity. 
Operation at other than the design point tends to decrease efficiency 


meee norten, markedly, the useful life of the device. 


iboats restriction on exhaust velocity is treated in the literature 
as a departure from the ideal conditions of constant exhaust power. 
The author treats this restriction from a slightly different viewpoint. 
In the succeeding section parametric equations for a constant-specific- 
impulse rocket will be derived which are analogous to equations (2-8) 
through (2-12). With these parametric equations one may compare 
directly the performance of a rocket with specified maximum power, 
under the two modes of propulsive control, i.e. constant-specific-im- 


pulse (CSI) or variable-specific-impulse (VSl) 
2.5 The Constant-Specific-Impulse Concept 


The propellant-optimal thrust program for a CSI rocket is shown 
subsequently to have two magnitudes, maximum thrust or zero thrust. 
In order to allow for maneuvering and for departures from the refer- 
Seer ajectory, réserve thrust capability is required. (For Vol rock- 
ets this problem does not occur since thrust is continuously variable. ) 
A method for controlling thrust magnitude which appears quite attrac- 
tive and is compatable with thruster design is to build up the propulsion 
ioitietrom a large number of small thruster nozzles of constant specific 
impulse. Proposals of this type appear in the literature but are not 


395 


analyzed in detail Control of thrust magnitude is obtained by con- 


trolling the number of thrusters in operation. 


If a nuclear energy source is used the power range is not continu- 
ously variable from zero to maximum power and may be limited in the 


number of times that stopping and starting is permitted However, the 
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energy source and the thruster units can be operated al their respecting 
optimums by placing an external high temperature resistor in the sys- 
tem which will radiate energy in excess of that required by the thrusters. 
Since the nuclear fuel constitutes only a very small fraction of vente 
weight, this procedure has negligible effect on the final mass of the 


vehicle. Figure 2-a shows the propulsive system schematic diagram. 


For a propulsion unit constructed in this manner, the total thrmeg 


of the vehicle is some fraction of the total thrust available. Thatis 


f = -nmec (2-19 
O 
f=-me (2-144) 
O O 
ine nm, (2-14b) 
nt, -nm_¢ 
aes Se eee (2-15) 
11 mM 
“nm ¢ 
p=np = ——— (2-16) 
= 2 


where n is the fraction of thrusters in use and te -m and p, are the 
maximum values of thrust, propellant flow rate and exhaust power, 


Pespectively. 


By eliminating the exhaust velocity, equations (2-14) through (2-16) 


may be manipulated into the form 





C1] eee (2-158 
m 2pom 


Integrating and expressing the result in terms of the mass ratio, one 
obtains 

Uf 

m f 

Vike= v= 2 f aS 
2D, O m 


a dt (2-Tay 





Observe that the factor ay 08 is the maximum acceleration possible 
for a given mass. Therefore the acceleration integral for a thrust- 


limited vehicle is proportional to the integral over the thrusting time of 
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the product of the maximum instantaneous acceleration and the instan- 
taneous acceleration. Designate this integral J* and observe that for 
a vehicle of specified power and initial mass, the mass ratio is a min- 


fai when J“ is a minimum 


i 
1c = 1 = 
J : s sane (2-19) 
a Se 
MR-1= = Jj* E20) 
p 


Equations (2-19) and (2-20) correspond to equations (2-9) and 
(2-12) respectively, and permit a direct comparison of the performance 
of a hypothetical vehicle with given power and given initial mass when 


controlled as a VSI vehicle and then as a CSI vehicle. 


The comparison is made by finding an optimal trajectory for a 
specified mission for each type of control and then comparing the mass 
ratios or the acceleration integrals for the two cases. This has been 
done for a simulated mission to Mars using techniques described in 


later chapters, and for the simple case of a transfer in field-free space. 
Zo) Comparison of CSI] and VSI Control 


From a qualitative point of view, finding an optimal trajectory for 
the VSI machine consists of finding that trajectory which minimizes 
propellant expenditure for a given constant exhaust power but with un- 
sometrained exhaust velocity. In the case of the CSI vehicle, the prob- 
lem is that of finding the trajectory which minimizes the propellant 
expenditure subject to a given maximum exhaust power and a given 


constant exhaust velocity. 


If both vehicles have the same maximum exhaust power and both 
control programs optimize propellant consumption then constant-specific-—- 
impulse is but a limiting case of variable-specific-impulse and 
cannot result in less propellant usage. Therefore the acceleration 
integral for a VSI transfer can be used as a reference for assessing 


the additional propellant cost of CSI transfers. 


ye 


In order to compare the two vehicles for a transfer in the srayitee 
tional field it is necessary to use machine computation. However,they 
may be compared analytically in field-free space. The results are 
shown to provide reasonable approximations for the values obtained in 


the gravitational field. 


Assume that in field-free space it is desired to traverse the dis- 
tance Lin the time T such that the rocket begins and ends at rest. This 
problem is solved in the literature for VSI spacecraft’ The solution 
is rederived in Appendix A and the soliition for the Csl spacecratma. 


also derived. 


The results for variable-~specific-impulse thrusting are: 


Z 
Scie ; 
tia ¢ Rae (oe 
a (2-22) 
ac 


where ao is the initial acceleration. The optimum acceleration pro- 


gram starts at a and decreases linearly with time such that a(T)=- ao: 
(see, Edcnrerz-c) 


ACCELERATION 





Fig. 2-c. Optimum acceleration for constant power rocket in field-free space. 


18 


For a normalized mass (i.e. mo is unity), equation (2-12) may be 


meritien in this case as: 








2 
Nie ee (2-23) 
pr 
feiime the parameter R. 
Re ch (2-24) 
p 
eT (2-25) 
3 
Therefore in this case 
2 
R = se (2-26) 
pT 


lation (2-25) is plotted in Fisure 2-d. Ris an excellent measure 
Sime cifficuliy of a mission in field-free space. For example: A 
Paeee iransier distance, a short transfer time and a small energy 
source would produce a large value of R- Such conditions indeed repre- 


sent a difficult transfer. 


li should be noted that the mass ratio is linear in R. Thus V5sl 


transfers characterized by a finite R require finite mass ratios. 


For a constant-specific-impulse vehicle of identical power to per- 
form this mission requires a specification of exhaust velocity or, more 
Semveniently, initial acceleration. In Figure 2-d, the mass ratio is 
plotted as a function of R for two values of initial acceleration. The 


first corresponds to ay = oa and is aSymptotic to R = 20. The second, 
az 


which corresponds to an optimum value of initial acceleration for CsI 


machines is exponential in R. 


iieveeliumuim initial acceleration for CsI vehicles is 


a = 


(In MR) a (222) 
lo t) O 
ect 


| (Opt) yoy 


I 


For small values of R such that the mass ratio is near unity, equation 


(2-27) may be simplified to 
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a, ee ao (2-28) 
(Optles, * = (opt) ysy 


The transcendental equations for this example are derived and solved 


in Appendix A. 


The conclusions that may be drawn from the field-free space 


example are: 


1) Control of thrust by varying specific impulse but with con- 
stant exhaust power is the propellant-optimal method of con- 


trol. 


2) For missions represented by small values of the parameter 
R, (less than 5) the propellant penalty for using an optimized 


constant specific impulse is less than 15%. 


poe For large values of R the propellant cost of constant-speciiic— 


impulse Is avery ldtce. Cven ton an Optimized Gon 


4) For sufficiently long transfer times the value of R may be 
decreased so that the propellant penalty of constant-specific- 
impulse is very small. (As will be subsequently shown, an 
optimized one way trip to Mars corresponds to an R of approx- 


imately 1). 


mois not readily apparent that the results derived for field-iree 
Space are directly applicable to transfers in the gravitational fields of 
the solar system. However, Melbourne and Sauer" have computed 
the Vol acceleration requirements for a number of interplanetary trans- 
fers and found the values to be in surprisingly good agreement with the 
Tield-free space analysis: Therefore, it seems valid to argue on the 
basis of field-free space when comparing CSI and VSI systems. The 


numerical results in this study support this conclusion. 


The parameter R, which includes the effects of both mission re- 
quirements and spacecraft power, was found to be more useful in the 
comparison of CSI and VSI systems than the acceleration integral alone. 
some significant facts which are evident when R is used as a parameter 


but which are missed when the acceleration integral alone is used are: 


ot 


1) increasing the spacecraft power for a particular mission 


reduces the propellant penalty of CSI operation, 


2) the optimum initial acceleration, and thus the optimum ex- 
haust velocity for CSI vehicles may be determined to reason- 
able accuracy from analysis of the VSI spacecraft using equa- 
ton (2-27. 


2.7 Maximizing the Payload 


In section 2.6, the problem of maximizing the mass ratio is con- 
sidered. In this section the energy source which maximizes the pay- 


load for a given mission and mode of operation is computed. 
Both CSI and VSI vehicles are characterized by the general equag 
ion: 
MR=— J+1 (22909 
p 
m 
= 


where the term B J may be written alternatively as ei 
p 


If the initial mass of the vehicle is considered to consist of only 
three parts: payload, propellant and power source, then the mass dis- 


tribution may be written as 


us eee | (2-30) 
m MR 
O 
m 
Soe ee B (2-3) 
m MR 
O 
pale 
where —— is the payload fraction. 
m 
O 


Substituting equation (2-29) into equation (2-31), an expression for 


payload fraction is obtained in terms ofa, 6, and J. 


/ 
-L . gf_1 - ) (2=02m 
m NawJoare 


Ze 


Differentiating equation (2-32) with respect to B and setting the deriva- 
tive equal to zero results in the optimal mass distribution for either 


@oor Mol space vehicles, The result is: 


m7 
(= | 4G ED (2-33) 





ae 
opt 
B= (a3) 12 -aJ (2-34) 
Opt 
Pri 1, 
p =(aJ) /* (2-35) 
70 
opt 
m 
where —P is the propellany iraction, 
m 
O 


From the alternative forms of equation (2-29) and equation (2-34) 


the optimum exhaust power for a space mission is obtained. 


Ls 
opt ors =) jay (2-36) 


Oi IW a 
O 


For the numerical work in this thesis an optimistic, but not un- 
reasonable, value of a was selected. In all subsequent computations, 
ana of 10 kilograms per kilowatt is assumed. From the digital com- 
puter studies of a 150 day Earth-Mars one way transfer, a value of J 
was obtained for the sizing relationships. This value of J was multi- 
plied by the factor 4 and rounded off to provide a more realistic value 
for around trip to Mars. The use of an approximation is Justified in 
that the purpose here is to provide a reasonable value of power for com- 
paring VSI and CSI operations without placing undue emphasis on opti- 
mizing the round trip mass distribution. That is, the author is more 
concerned with optimal guidance of a space vehicle of given power and 
mode of operation than in determining whether the spacecraft is the 
best vehicle for the mission. Further, the preceding method of maxi- 
mizing payload, although widely used in low-thrust-studies, does not 
adequately treat a round trip mission comprised of several phases such 


as escape from Earth, transfer, capture at the target, a wailing period, 
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and similar return phases. Each of the phases may impose particular 
requirements which are not compatible with the other phases. For 
example, it may be impractical to use 100 or more days to effect a 
spiral escape even though the optimized solution requires it. Conse- 
quently the derivation in this section must be considered as only an 
approximate method of maximizing payload for realistic round trip 


mission planning. 


The value of J selected was lone The units are (astronomical 
units)- per (day)”. Inserting these values of a and J (in compatable 
units) into the equation (2-36), a value for p/m, was computed for use 


in the numerical work. 


Po 


0. 0242 kw/kg 


m4 (27S) 
3 


6/6988 X10 “(Ac \-fiday) 


Using the value obtained in equation (2-37), values of MR were 
extracted from computer results and compared with the field-free 
space values for both CSI and VSI engine control. The results of this 
comparison are plotted on the curves of Pisure 2-d, ) sereeimecnia 
good for the cases tested and warrants further investigation for dif- 


ferent values of the parameter R. 


Values of initial acceleration and specific impulse which are char- 


acteristic of the values obtained in the numerical studies are 


+ 


a /8, eae (2-38) 


J 4000 seconds (2-35 
Sp 
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Clare TE i ii 
LINEAR GUIDANCE 


Seieouimimary of Chapter II 


moe suidance problem for continuous-ithrust vehicles is formulated 
in this chapter using perturbation techniques and the theory of optimum 
control. A reference trajectory which satisfies mission requirements 
is assumed to exist and its associated control program is known. ‘The 
eemsequences of allowing the reference irajectory to be an optimal tra- 
iee@rery are examined. Linearized equations of state are solved to- 
gether with their adjoint set to produce a state transition equation 
applicable for small perturbations around the reference trajectory. 
The state transition equation is then shown to be suitable as a guidance 
equation. Solutions for the guidance equation are derived for fixed- 
Beeeot-arrival (FTA) guidance. Variable-time-oi-arrival (VIA) suid- 


elec is solved for certain restricted cases. 
3.2 General Remarks 


In terms of the discussion in Chapter II, the ensuing derivation of 
ita: ouidance may be characterized as a method of minimizing the 
eects avion integral, J, between the present position and the target 
position. For this purpose it is convenient to abandon temporarily the 
engineering aspects of low-thrust transfers and to consider the prob- 
lem in terms of the calculus of variations and the theory of optimum 


Control, 


The guidance schemes suggested for interplanetary vehicles fol- 


lowing ballistic trajectories usually require that midcourse corrections 


Siare Veansitlon Cquation, An equation which expresses the state of the 
vehicle at any given time in terms of 1) the state at any other time and 
2) the control existing between the two times. 


guidance equation: An equation, possessing a solution which satisfies 
the mission requirements for the spacecraft. The solution of a guid- 
ance equation is oftentermed the "control law" or "control program" 
for the physical guidance system. 
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minimize, subject to one or more constraints, the position and veloc- 


ity variations at the destination point. 


The characteristics of a continuous-thrust vehicle permit a curds 
ance scheme which will null the position and velocity variations at the 
destination, provided guidance is begun sufficiently early in the trans- 
fer. For all practical vehicles, if guidance is not initiated early imma. 
transfer the low-thrust vehicle will arrive at some point, beyond which 
insufficient power is available to complete the mission. The vehicle 
is then said to be ''in extremis'. This boundary may be visualized as 
a conical surface surrounding the reference trajectory with its "apers 
at the target point. The ''distance" in phase space from the trajectory 
to the surface depends upon the excess power allowed for guidance. 
The problem of determining the optimum power and propellant allow- 
ance for guidance is not studied in this report; however, the techniques 
of this caapter are readily applicable to such investigations and are 
recommended as.a tool for future work. In this chapter it is assumed 
that the reserve power is adequate for satisfying the solution to the 


guidance equation. 
3.3 Formulation of the Problem 


The basis for the subsequent analysis is that a vehicle is in transit 
between two planets. The control program in use corresponds to some 
known reference trajectory. Measurement of vehicle state indieares 


that the vehicle is off the reference trajectory by some small amount. 


The first problem of interest asks the question: ‘'What is the 
effect of the state variation at the time of the measurement on thew 
hicle state at any future time?'' The next obvious question is: ‘What 
is the new control program which will cause the vehicle to satisfy 
mission requirements?'' The third is: "If there is more than one con- 
trol program available, what criterion should be used for selecting one 


of them?" 


The third question may be disposed of immediately. The thesis 
is concerned with propellant-optimal guidance. Therefore the choice 
between controls is on the basis of minimum propellant consumption. 


Prom Chapter il, aeriterion which satisties 1his objective as the 
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Peeeclcration integraly There may be other criteria which are equally 
as good, however J was chosen because it is easy to manipulate and is 


Eecily amierpreted in the physical problem. 


miessecond question requires that the mission be defined. For 
this purpose it is sufficient to require that the final position and veloc- 


ity of the vehicle have certain well defined values. 


Depending upon the particular mission, the above guidance object- 
ives may be satisfied at the preplanned arrival time (FTA guidance) 
or at some time different from the planned time but consistent with the 
Mmipssion (VITA guidance), The FTA problem is readily solved by the 
linear methods of this chapter. The VTA propellant-optimal problem 
Peeconsiderably more difficult. A method of obtaining restricted solu- 


tions by linear methods is presented and discussed 
o.4 selection of State Variables 


The differential equations of motion for celestial bodies admit six 
Somstants of integration. The selection of the six quantities to repre- 
Pem@e this motion is to some extent arbitrary. To be consistent with 
ihe definition of mission requirements in this investigation, and be- 
Beuec they are convenient, the three components of position and the 
three components of velocity are chosen. The convenience is due to 
their facility for visualization and their relative facility for physical 


measurement. 


Since the conservation of propellant is important, an additional 
differential equation describing the mass change due to propellant flow 
must be included for certain cases. Thus seven independent, but 
coupled, variables are sufficient to describe the state of the vehicle 
at any time. The state will be subsequently written as a vector which 


has components of position, velocity and mass. 


Es 


Se ai, 
m 
The second order differential equations of motion will be rewritten as 


first order equations to conform to the above definition of state. 


2k 


For this thesis, only central force fields are considered. It is 
subsequently shown that extension to more general gravitational fields 
does not invalidate the theory but does introduce computational com- 
plexities which serve only to obscure the physical concept if introduced 


Pre tia pumel.. 


With this introduction the differential equations of state are now 


presented. A nonrotating frame with its origin at the central body is 








assumed. 
S = g (s, a) (3-2) 
where me 
[ 
gid Krprtap (3-3) 
— r : 
6 ti 


and where a is the thrust acceleration, subsequently called the Comma 
vector, W is the gravitation constant for the central body, and g_,(m. a4) 
is the mass rate equation for the particular vehicle under considera- 


HO Leon @ ta prcrr al 





ae 
g_(m,a) = - ——— (VSI) (3-4) 
2p 
gta) = - = (CSI) (3-5) 
ce 


Throughout the report the variable m is a normalized mass, that is 
m= 1. This obviates the necessity of selecting a value for totais 
hicle mass. The power and the propellant flow are likewise normalized 
variables. Thatis 


_ Exhaust power 
2 eee eee (3 =o) 
iitloa! macs 


ee propellant flow rate (32m 
initial mass 
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peo linearization of the Equation of State 


The assumptions of a known reference trajectory, known control 
pHoOerammand Of small deviations from the reference trajectory permit 
the application of linear perturbation theory. Thus 
og 


6s +— 6a (3-6) 
oe 


From equation (3-8) the matrices A and B are defined. 














? on ae 
A as “G ©, © Cec 
0 9 
a of of %m 
om 
ro = 
dg 3 
Be— =} 1, (sao) 
Gal 9 
ae 
where 9 / i 
—— for (3-11) 
or ie ? 
G = 4 ee eae eee) 
a fee 
oe ae a 
a (VSI) (Bas) 
om p 
98m _ a 
Se see: (CS) (sai 
om e 
dg Zo 
eee voi) (3-15) 
da p 
ne 
ae iy. 
= EST) (3- 15a) 
0a GC a 
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The matrix Gis a symmetric matrix of gravitational gradients. The 
extension to a general gravitational field does not change this property>” 
However, numerical integration in such a field requires an ephemeris 
for the disturbing bodies. This is a computational problem and does 


HOtameci ine tCOry, 


Equations (3-8) through (3-15a)treat the variations in state velocity 
as functions of the variations in state and variations in the control vec- 
tor, a. It is often desirable to study the effects of variations ingame 
formance parameters such as power, exhaust velocity, propellant flow, 
etc. For this purpose one may rewrite the vector a and the propeliam 
flow equation in terms of the appropriate parameters. New matrices 
A and B will be formed for this purpose. These cases are deriveqdmm 
Appendix C since they do not contribute to the discussion of guidance. 
However, occasionally in the guidance problem it is desirable to work 
with the thrust, f, as the control instead of a. The theory subsequently 
presented is applicable to this formulation also, Nor the presen 


however, the form 
6s = Ads + Bba (3-16) 


is considered to be the fundamental formulation for the guidance prob- 


lem. 


With the perturbed equations of state, as defined by (3-16) , and 
using the adjoint method, it is possible to derive a state transition equa- 


OMe 
3.6 Method of Adjoints 


One description for the method of adjoints is obtained by consider- 


ine a set of equations related to4%3- 16) by she mari cumctien 


=- AA (3-17) 


Equation (3-17), with arbitrary boundary conditions, is said to be ad- 
joint to (3-16) and the elements of A are the adjoint variables. If equa- 
tion (3-16) is premultiplied by A , (3-17) post-multiplied by 6s and the 
resulting equations are added, then one obtains 


4 ( Abs) = ABéa (3-18) 


a0 


Pwiemed Siven reierence trajectory, i.e a = OQ, equation (3-18) has the 


solution 


Abs = constant vector (3219) 


feo 1) is integrated along the trajectory, subject to suitable boundary 
conditions, then A. is a known function of time and the state variation 
at any time may be determined from the state variation at any other 


ime Provided the reference control program is used. 


The adjoint method may be described as the process of introducing 
a set of known auxiliary variables, A. , which satisfy (3-17) and which 
transform the state variation at a given time into the state variation at 


mipeouner time, provided a forcing function does not exist. 


A scalar approach to adjoint equations and other examples of their 


application are presented in Appendix B. 


A useful system of equations closely related to the adjoint set is a 
set often called the fundamental set or fundamental solution. It satis- 


fies the relation 


b= Aad (220) 


Combining (3-20) and (3-17) by the appropriate pre-and post-multiplica- 


tion, adding and integrating, yields the first integral 
AD = constant matrix ead 


If the boundary conditions are chosen such that 


Ait) = 1 (3-22) 
O(0) = 1 (3-23) 
then eee Oe H(t.) (3-24) 


This relationship is useful in later discussions. 


3.7 The State Transition Equation 


If equation (3-18) is integrated between the times ty and to, jie 
result is 
t2 
oe ssi J ABéa at (3-25) 
1 
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Then, provided A 9 is nonsingular 


2 
_ acl -1 ; 
OS 5 wen Coen - AB &a dt (3-26) 
1 


Aeprool 1s presenled im A peemelent] cman yee is not singular. 


Tice reduce a N, is called the state transition matrix and is 
denoted by aa where i and j] represent any two times. Equation (3-26) 
is called the state transition equation. It transforms the state varia- 
tion at ty and the control perturbation between ty and to into the state 
variation at to. 


There are three useful applications of the state transition equation: 


1) It may be used as a tool for studying the sensitivity of the 
trajectory to perturbations in launch conditions and to anoma- 
lies I engine periat mance: Friedlander!® has performed 
some investigations in this area. This application is not 
pursued further in this report except for derivation of apelims 


cable formulations in Appendix C. 


2) If engine performance is Measured by an accelerometer ion 
other device which can be related to the state equations, then 
the state transition equation may be used in the navigation of 
a spacecraft by relating the measured engine performance to 
the state of the vehicle. This application is discussed in 
Chapter ly. 


3) The most important application of equation (3-26) is its use 


in guidance. If, in (3-26) , t, is considered as the final time, 


then the state and control ee occuring along the trajectome 
may be related to the final state variation. Since the mission 
objective may be defined by the final position and velocity vee= 
tors, (3-26) meets the requirements for a puidance equation 
provided that, for any state variation Os | it is possible to fina 
a solution a + 6a which will cause the vehicle to match the 
physicalboundary conditions at the terminal point. For this 
application the terminology ''guidance equation" is preferable 


to ''state transition equation" and will be used to define the 


oe 


equation when t, is the final time and the boundary conditions 


on the adjoint tere are chosen as A_(t¢) = I. Thus the 
guidance equation 1s 
ne 
bs, = A, esyt J A.B 6a dt eeu) 


where the initial time t is any time between the launch time 


and final time, at which a new control program is desired. 


itie integrand in (3-26) and (3-27) is observed to be quite simple 


when the matrices are partitioned and expanded. For the control vec- 


ren 2 
cee ee 
N=) Ao1 Azz © eens) 
alt ae 
oF O° Ag, 
Using (3-10) and expanding 
N12 
> oe 
A.B TOD ( ) 
Jean 





where /\ 33 is a scalar and the remaining J's are three by three 


mieirices, Jhe vectors are three component null vectors. 
3.8 Solution of the FTA Guidance Problem 


The arguments presented in this section are equally valid whether 
the desired solution is the acceleration program a+ 6a or a thrust pro- 
gram f+ of. A discussion of the requirements for and the consequences 
of interchanging control vectors is presented in Appendix C. Since the 
propellant expenditure is so easily written in terms of the acceleration 
integral, especially for VSI vehicles, it is computationally convenient 
to work with a as the control. If ais used in solving the guidance prob- 
leis vic set OF adjoint functions, A, may De Seduced limimediately Irom 
a seven by seven matrix to a six by six matrix since the equations of 


motion, (3-3), do not contain mass explicitly and the mission" is to 
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satisfy terminal conditions on the equations of motion. However, when 
working with the CSI vehicle it is usually more convenient to use the 
thrust, f. Acceleration is used for both vehicles in this report in Gmdam 
to be consistent. Since thrust and acceleration are relaved tire 
mass, the equations of motion using thrust as a control must contain 
the state variable, m. For this reason the adjoint set cannot be re- 
duced a priori to a six by six matrix. In recognition of this fact and 
since both control vectors are useful, the A matrix is always con- 
sidered to be a seven by seven matrix; with the understanding that if 
ais the control, certain elements are identically zero. This doesmiai 
present any problem in manipulation and need not be considered fur- 
ther except when A. and B are to be evaluated numerically. In the 
subsequent discussion, however, a reduced adjoint set is introduced 
and is denoted by J_*, This reduction is not associated with the dif- 
ference between acceleration and thrust formulations. Both a A ma- 


trix anda A. * matrix exist for each control vector. 


From equation (3-27) it is clear that if a state variation exiSitcees 
time t, then in the absence of control changes, a variation in final 


state will occur which is given by 


5 = z 
s,= A, 6s, (3-28) 


The guidance criterion is that the final position and velocity must satisfy 
mission requirements. Presumably the reference trajectory satisfies 
these requirements. Thus the position and velocity variations fronmmmae 
reference must be zero at the final time 1 FP PAs cuidance tosses 


Stated mathematically, it is required that 


Osta | O (3-29) 
5m 


where Om is a small but unknown variation in final mass. Whatever 
the value of Om, it is of no immediate concern. To reflect this itis 
convenient to drop the 6m, and define the ''miss' at the target which 


arises from the state variation at time t as 


ow / N, Os | (3-30) 
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where € is a six component vector of position and velocity variations 
and /L*is a six by seven matrix obtained by deleting the seventh row 
of A.. Clearly 6m 


known, from 


, can be obtained, once the corrected control is 
t 
_ ly ae 7 
fm, - Ao, b8,+ J NN, B 6a at (3-31) 


where UX. ~ denotes Phere ltcdeoe VERE Lon Or iL Mewad|Olmie sci: 


since the state error 6s, is assumed to be small, so that linear 


theory is valid, it follows hat Sesidal wCOngec lon tenniic melereiocEGe nl 
ieeimorogram will be sufficient to null the miss vector in the remaining 
ment time. Lhe foregoing is true provided the vehicle has sufficient 
thrust. Thatis, it is not 'in extremis''. Assume it is not. Then the 


miss may be reduced at a rate such that 
£= A*Béa (3-32) 


or by integrating, such that 


O-é& = f A* Béadt (3-33) 
That the control 6a is not unique, except for a gy which requires 
maximum thrust continuously, is evidenced by the fact that for ballistic 
guidance, ideally, only two corrections are needed to null position and 
melocity error; a midcourse correction to correct position and a termi- 
Maleeorrection for velocity. Thus for continuous thrusting vehicles an 
fie y Of Solutions exists. The criterion tor selecting a unique con- 
trol has already been given; the acceleration integral J or J* must be 
Minimum. Mathematically, the problem for VSI vehicles may be 


BuaLved: seuss required that 
t 
i 
o| ee. a) AB fa dt (254) 
O | t 
and that J is a minimum, where 


t 
j= - (a + ba)! fay Oa) 


= dt (3-35) 
t 


pee 
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; O 
For convenience define the optimal contro. 


a a (3-36) 


For a VSI vehicle with no constraints on a> the set (3-34) and (3-35) 
can be solved by a direct application of the calculus of variations. The 
vector space around the reference trajectory is explicitly assumedsae 
be ilat toairst order, Consequently /\. and B are invariant between 


neighboring trajectories. Thus 


Ur of o te 
§ (J =F see Ee l + ‘ vas: Bia® - a) sf) = 0 
t 2 2 a t "Se 


J 


where 7 is a vector of constant Lagrange multipliers and the subscript 


(3-3 


on om is dropped. Expanding (3-37) one obtains 
et - ib . 
a =-a7 A*B (3-38) 


The vector 7 is eliminated wali melo er eqnanlen (> >55))r 


. Ur 
| | = &€ - A* B Be Axl i a dt 3-39 
g f A ( ) 
t 


jof 


a 
-7=M + (n-€) (3-40) 
where 
Mi 
m= f AxBB! Ax! at (3-41) 
f 
4 
ie ene eee (3-42) 
Thus 
a° = B Ax’ Mm? (ny -é} (3-43) 


The solution (3-43) is valid provided the reference trajectory 1s Sutme 
ciently close to an optimal that 6a is small. It is unique provided M is 
Not Singular. A proof is presented in Appendix C for the Cxistenceges 
ie | 

A physical interpretation of the quantities M and 7, defined by 
(3~41) and (3-42), is valuable in understanding the solution. First ob- 


serve that € has two interpretations. From (3-30), € is the miss 
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arising from a state variation, OS 4, MWiciiiiewtCtreMeCemcOuLrol tS ViSed, 
From (3-33), - € may be interpreted as the miss that results from a 
control variation, 6a, occuring after time t but with os, =U rota c ibe wel hy 
mario -4+2), - 7 is interpreted as the particular miss that arises if 
6a = -a; that is, coast is initiated. If the miss, €, due to a state varia- 
tion and the miss, - 7, due to initiating coast are such that Cae) 
the optimal control is the null vector. Although such a solution will 
generally violate linearity assumptions, it will indeed result in a mini- 


mum J, namely zero. 


If the product A.* B is interpreted as the sensitivity of the final 
miss to a unit impulse in each component of the control at time t, then 
M may be considered as a weighted total sensitivity of tne miss toa 
unit control applied continuously between t and te. In a physical sense 
meewprocedure computes the vector sum of the miss using the reference 
control vector and the miss using a null control vector; then selects 
the control at each point according to the sensitivity of the miss at that 


Point. 


It is interesting to note the solution which results if, instead of J, 


Piesiminimization criterion is 


She eee (3-44) 
t E 
Proceeding as in equations (3-37) through (3-41) but solving for 6a, the 


result is 


ao=atbaz-a-Bi A*' mM’: (3-45) 


Comparing (3-45) with (3-43) it is observed that for a = Be Axt M +n 


Peewresults are the Same. ‘The preceding expression for a satisfies 
(3-42) therefore for small variations around the reference trajectory 
the two methods result in the same control. This is only true for VSI 


eajecitories, however. 


Before proceeding to consider CSI solutions it is desirable to 
Serine ne COUSeQucHees Ol Ine relecrence trajcerory beim Optimal 


for the initial launch point. If the reference is initially optimal but a 


of 


state variation occurs such that a miss, &, will result from use of the 
reference control program, the reference control ceases to be an 
optimum because the boundary conditions are not satisfied. Also the 
particular formulation used for the problem does not result in a singu- 
lar M matrix as is often the case. Therefore, in this study, the omg 
consequence of using an optimal reference is the assurance that 6a 


will not violate linearity assumptions. 


A new approach is now presented which can be used for both VSI 
and CSI control. To illustrate this method the VSI problem is solved 


again but with a constraint. 


Assume that thrust is not to exceed an amount tS: Pie rete. 
Mirust acceleration caniict exceed f./mm.. To be consistent with the 
normalized mass variable, m, used in this report denote hy as the 


initial acceleration limit ao The constraint may now be written as 


a< (3-46) 


Blo” 


Following Reine and using Pontryagin's principle (Appendix D), 


form the scalar Hamiltonian, H 


a 
O O 
a a 


oe ee (3-4 
Z 

The theory states that if the cost,|(a°)?/ 2Mtis to be a minimum, thenieg 
each point along the trajectory H must be a minimum, where v is an 
unknown vector, often called the costate, which satisfies the adjoint 
relationship, and & is the state velocity. The linearization in this chap- 
ter permits considerable simplification. Using (3-32), E is interpre 
as a variable representing the velocity of the final state. The varialgie 
fy is the initial condition for é and is a function of the lower limit of 
integration, t. 54 is evaluated from the vehicle state, OS 4, 
Because @ is a variable in state space at the final point, vy is a constant 


using (3-30). 


vector. In particular it is the final value of the general time varying 


costate vector. 


Ineretore, Using (3-36) 
Gao Se 8 (aa) (3-48) 
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Further 


ae 
O° Oo 
a 6a T O 
ee wee) (3-49) 
Z 
ale 
a Bias 
H=-v? A*Bat=— +y7 pepe? (8-50) 
Zz 


The first term in (3-50) is a constant independent of ae, thus His a 


MmiiMiauin when the sum of the last two terms, H,, is a minimum. 


ie 
The control a is to be determined such that H, el 2 igahlianbagihaer 
ar O 
wa a a O 
minimize H, = ——_—— +v /A*Ba (s-51)) 
2 
26 
oj cr Le a< — (3-46) 
=" «iy 
now let vi A* B=q! (3-52) 


pumee gq and ae are both three component vectors, one may be obtained 
from the other by a scalar multiplication and a rotation. Thus 


ee eel (3-53) 


where y is a positive scalar and C a coordinate rotation. Inserting 
feo) into (3-51) 
ee 


H, -“_— tyq' Cq (3-54) 


Herany value of y, H, is minimum if qi Cq has maximum magnitude 


i 
and is negative. But 


la’ cal < lal leq! (3-55) 


where equality holds for C =I. It is apparent that Hy is a minimum 


only if C = - land (y 7/2 - y) is a minimum. 


é 2 
H, = (2 q , (3-56) 
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Thus when constraint (3-46) is applied, Hy iS a minimum if and only if 


O oe ea 
a =_-— —— Q > —— (3-57) 
Eimer) 1a 
O 26 
a° = -q q< (3-58) 
m 
a 
y= min 1; —— (3-99) 
Inq 
oie 
ae eye eee ay (3-60) 


Comparing this solution with the variational solution, equation (3-38), 
one observes that except for the constraint y the solutions are identical 
and V=a. y may be interpreted as a switching function which Caria 
the thrust restriction (3-46). 


Having gained confidence in the Hamiltonian, now apply the method 


to the CSI transfer for which the classical approach is degenerate. ‘The 








appropriate cost is t¢ 20 
pa) eerie (2 =1ge 
iL 21m 
Minimize O 
aia T - 
ne +v~ A* Bia -a) (3-61) 
2m 7 
a 
subject to a< ae (3-46) 
mm 


Again vi A bra Ts AecOnsianeimerciore 


a 


minimize H, ae a + q? an (3025) 
Zataral 
Where q is defined as before. 
The previous argument holds, to yield 
a =-74q (3-63) 
and 
ate 
Hey (52 - a) 4 (3-64) 
l Zia 
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In this case H, is a minimum if and only if 


li 
ay =-y Bees v (3-65) 
wiiere 
26 
ae q —— 
au eee 
45 are 
Ve Reet SA eee 
mq Zin 


For the case q = a ,/2m, yis actually indeterminate. However this is 
not an important consideration for this thesis since equality holds only 


Poesiniinitesimal time periods. 


Equations (3-65) and (3-66) yield the 'bang-bang" solution char- 
acteristic of optimal trajectories for which the cost is a linear function 
of the control. Again, the solution is strictly valid only for perturba- 
tions around a reference trajectory. In order to complete the CSI solu- 
Momit is necessary to evaluate the constant vector v. The evaluation 
fmemore difficult in the CSI case than for VSI because the control pro- 
gram is discontinuous. Vis evaluated by first applying the boundary 
conditions (3-34) and using (3-36). 

uF 
O=&+ {| A* B(a® - a) dt (3-67) 
IU 


Assume that ae eMcsa Citlemeonihy beeause pei Respective values on 7 
pide -are different, The assumption merely implies that /.* Bis the 
same for neighboring trajectories. With this assumption rewrite 
(3-65) as 

ao = - 7° Bax! y° (3-68) 
ig@jne corrected control program, and let 

BL aay Baie v (3-69) 


MemGcsenb une relerence control procram, Phen (3-6/7) becomes 


t 
f 
O=£- f A* BB AK! (72 v°-yv) at (3-70) 
t 
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which may be rewritten as 
Uf 
SS f as Be Axl A(y v ) at (3-71) 
t 


The increment 4 (y vy) represents the difference am (he control soem 


gram between the reference trajectory and the corrected trajecromm 


since y is a function of v, a solution may be obtained by formally 


Giilerentiating and solving (j=72) or oa, 


i 

is 

-|29. f AxBpBt Axtyvat | av (3-72) 
Ov t 

The functiony is discontinuous at switch points, which occur at times 

ty in the interval te =t. Conseqtently the intesral{s=/Z) mice 

separated into regions of coasting and thrusting and Leibniz's rule 


used for the differentiation. A term of the form 


eee ot, 
+ A* BB A* y v oe Ov results for each switch point, where 
OV 


plus is for switch off and minus is for switch on. 


The term dy /8v is evaluated in the continuous regions from the 


definitions (3-52) and (3-66). The term ot, /9 vis evaluated by con- 





sidering 
q=q (t,v) (3-73) 
Proms 06) 
a 
q(t,.¥) =— (3-74) 
2m 


Differentiating (3-74) with respect tov, one obtains 





Ot 
CGE Hee eG (3-7 
ot OV OV 
0 
Ot Ov 
Me ss (3-76) 
a oq 
us Ot 


The derivatives of y and q are evaluated in Chapter VI when it becomes 


necessary to specify the explicit form for computation. 


Solving for Av, one obtains 


Av=mM* te cone 
where 
“f T. OT 9 TT ane 
M* = [AX BB’ A+ ly + or ee ea Bas e/a ay , a 
c Op k av 
a (sere) 
The constant vector vo is obtained from 
vo =p+t+Ap Ga79) 
where 
vy=-M!y (3-80) 
tf 
M= fy A*BB A*? at (3-81) 
7 
enc 
Ms 
n= { A® Ba dt (3-82) 
t 


mae values for y° are obtained from equation (3-66) by using the com- 


puted De ied. 


Witn all quantities defined in the preceding equations, the corrected 


control program is given by 


2° = 7° BT AST (wety- Meg ) (3-83) 
3.9 Application of the Guidance Theory 


imporder to use the theory of section 3.6 for the guidance of space 
vehicles it is necessary to compute the quantities in equation (3-43) 
for VSI guidance or in equation (3-83) for CSI guidance. Since (3-43) 


may be regarded as a special case of (3-83), only the latter is discussed. 


[ietem Clete VechlnewCeClorvais KMOWieiel tie Clements Of tae ima 


trices B, A., M, and M* and the components Of (ae vector 7) can be 
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computed for any point along the trajectory Whether these quantities 
are precomputed and stored prior to launch or are computed on board 
as needed will be determined by the state of computer technology when 
the vehicle is designed. This problem is not relevant to the presem 

discussion. In either case, the first step in computing the corrected 

control is to establish that a state variation, OS |. exists. Methods for 
processing measurement data for this purpose are discussed in Chap- 


ter 1V. Then the fimal stateverror ie cerermincd (rom. oir 


B= At 


t 8s | (3-30) 


iie vector v° is computed from the known quantities and the final state 


error USING (3-1 i) neue 32. 





vo=mM*x'é -M ty (3-84) 
Fronts -92) ands = 06) q° and y° may be computed. 
g° = [BtAxt 7? | (3-52) 
O O oO 
a! Sa | 
os (3-66) 
a a 
OG O O 
ye ee { 
mq Zin) 
Then from (3-68) 
ao a ee Bi Ax? vy? (3-68) 


This result may be programmed into the vehicle control system to 


implement FTA guidance. 
3.10 Discussion of the VTA Guidance Problem 


In Chapter iJ it is shown that in field-free Space, the minimum 
acceleration integral for VSI vehicles is J = ee A VTA guidance 
scheme which minimizes J for a given L in field-free space wiliinem as 
fore select an infinite transfer time unless constrained. If the planets 
were all in coplanar, circular orbits and the reference trajectory lay 
in this plane, a similar result would be obtained in the solar system. 
Because of the inclination and ellipticity of planetary orbits, local min- 


ima of J will occur which depend upon planet and spacecrait 
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orientation. A complete relaxation of the time constraint, te, in the 
theory of sections 3.8 and 3.9 would require that the new arrival time 
correspond at least to a local minimum of J. Therefore, unless te is 
approximately equal to this propellant-optimal arrival time, linear 


theory cannot be used to obtain propellant-optimal VTA guidance. 


bevever, 11 1s possible io restrict the change in arrival time, use 
linear theory and obtain useful solutions. A case where this is impor- 
tant occurs for a vehicle "in extremis'. That is, insufficient power 
(or thrust) is available to null & at time te. The result obtained is 
dependent upon the method of restricting At and upon any simplifying 
assumptions. Consider the following example which is based upon the 


assumptions: 
1) The target point relative to the planet center is unchanged. 
tye he final velocity relative to the target point is unchanged. 


Treo inist acceleration a = 2 (t,) is constant Over the intervals 


i 
The motion and position of the target point at the time te Te eetee 
re (tp + At) = Fm (ty) ee At eS Ste) 
ae ae ee ee a At (3-386) 


where ae ida pete lle Diehectany IMeon anon tieidesateqcelaiive 


R 
velocity, respectively, and gr is the solar gravity at the planetary 


eects , 


For the mission to be accomplished, the vehicle position and veloc- 


apes! equalr andy, attimet,+ Ai. Phe vehicle position and 


ay aly ie 
wellocity in terms of the reference trajectory are 


Z 
r(t,+ At)=r (t-) + or (tp) + v (tp) At +a, =) + ov (tp) At 


(3-87) 
v (tp, + At) = v (tp) + dv (tp) +a, (At) +g, At (3-88) 


polving equations (3-85) through (3-88) for or ¢ and OV ¢ one obtains 
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Oa Ve nee en a 


i R aes (ST ee 


ov 


Ve=7 a At (3-J0pe 


aor 


For the linear approximation assume 67, /2 Say Equations (3-89) 


R: 
and (3-90) may be inserted in the guidance equation. 


(y ( ‘f 
-£-R\at = + ff A*B ba at (3-91) 
ar | = - 


Define a new miss vector &* such that 
] —- 


| ae oie Vv 
= { a 4 : At (3-92) 
5, Sy —f 


Then by selecting Ot such that the uncorrected final position errom 
|e 


time te + At and with the least departure from the actual irajecions 


| 
5 


fp 





, 1S a Minimum, the vehicle will attain the terminus but at the 


The result is 


pee (3-93) 


Equation (3-93) may be inserted into (3-92) and the optimal control 


program found using the PTA procedures of S€cilons)s,o andao.: 


It the target point is the planetary sphere oi influence and the des 
sired relative velocity is zero, obviously the above solution is invalid. 
Other assumptions may be used to treat such cases. The assumptions 
to be used and the criterion for restricting At may be changed te cum 
the purposes of the investigator. In general, the propellant-optimal 
VTA problem is not readily treated by linear methods unless additional 
constraints are used. Exploring the numerous possibilities that arise 


will be left for future investigations. 
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Cre ee iy 


Doin Om THB STATE VECTOR 


Pelee introduction 


In this chapter a method of determining a "best estimate’ of the 
vehicle state vector is presented. The term ''best estimate" is used 
to designate the estimate for the state vector computed by processing 
redundant measurement data through a statistically optimum filter. 
Potter and Stern’? have shown that all of the three commonly used 
methods of processing redundant data. namely. maximum likelihood, 
minimum error ellipsoid, and minimizing a characteristic scalar 
parameter all result in the same filter. They have shown further that 
for each unbiased optimum filter, there exists an associated biased 
optimum filter which produces a smaller error ellipsoid for the esti- 
mate . 

Since all of the methods result in the same filter, in this thesis 
the method which presents the least mathematical complexity is used. 
The method involves a slight variation of the minimum error ellipsoid 


ieecnnique . 
4.2 General Remarks 


To establish a corrective thrust program ii is desirable to have 
accurate knowledge of the state vector at the time the program is to be 
initiated. It is clear that given perfect knowledge of some prior state 
weeror and of vehicle performance the prediction problem reduces sim- 
ply to application of the state transition equation between the time for 
which the state is known and the time for which the prediction is de- 
sired. It is equally clear that perfect knowledge of a precess seldom 


exists. Thus the problem becomes that of using existing information 


An unbiased filter will produce a true value for the state vector if 
all measurements are free of error. A biased filter is biased in favor 


of the a priori or prelaunch expectation of the second moment of the 
state vector probability density. 
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in such a way that its inaccuracies have minimum effect on the guidance 


decisions. 


For ballistic vehicles the determination of position at two points 
along the trajectory suffices to specify the entire trajectory. Several 
investigators including Battin’” and Stern?! have studied the eifeers 
of using additional fixes and of selecting the geometry of individual 
celestial measurements to reduce the error in the trajectory deterimames 
tion. The techniques developed in those studies may be extendedwam 
some cases without change, to the present work. The primary differ- 
ence between the ballistic and continuous-thrust vehicles. is that the 
future portion of the actual continuous-thrust trajectory can never be 
completely determined on the basis of its history alone. This fact is 
due to the possible occurrence of random changes in the thrust vecigim 
Pnother diiterence whieh is more Sasily treated iastme ie sig@ wees 
dependence oi the trajeciory upon vehicle mass. This dependencemmag 
be handled by measuring propellant state at discrete intervals as well | 


as making celestial measurements. 


Thus for the continuous-thrust spacecraft, state determination from 


on-board measurements may be separated into two related problems: 


1. determining the state history by suitable filtering of all meas- 


urements that have been made, and 
2. predicting fulure values Or icrsrate, vector 


The second of these is dependent upon but is not uniquely determined 
Dy tne first, 


The use of celestial sightings, as in the case of ballistic transfem 
is sufficient to determine the spacecraft trajectory. However such 
measurements are not sufficient to specify the entire state history nor 
the variational history of the control vector which contributed to the 
state change. Because of this, prediction of the state vector solely 
from periodic celestial sightings and propellant measurement will be 
subject to larger uncertainties than if the prediction includes meas- 


urements of the control vector as well. 
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4.3 Unbiased Estimate of Present State 


The term "present" state will be used to denote the state vector 
computed without regard to a computational time lag, thus, it is the 


state at the time of the most recent measurements. 


The measurements which will be processed for the estimate are 
1) periodic determination of position, 2) propellant state, and 3) con- 
tinuous measurement of engine performance. It appears quite feasible 
to derive a method of including raw celestial observations in the compu- 
tation process without explicitly deriving the position vector. How- 
ever, this would contribute nothing to the present argument and will 
Memeit aS a Subject for future study. Consequenily, the computed com- 
ponents of position variation will be considered as a ''measurement . 
meee velocity is impractical to measure directly, except near planets, 


the velocity will be derived from successive position measurements. 


The time of the present state estimate will be chosen as occuring 
between celestial ''fixes''. If the state at the time of a fix is desired, 
such an estimate will be a special case of the more general problem. 
This approach is justified on the basis that celestial observations may 
be separated by several days but current engine data is always avail- 
eae. 


In the following development measured quantities will be denoted 


by the "tilde" (~) and the estimated quantity by a carat (A). 


The relation between the measurement quantities and the state may 


be written as: 
Nos = oq (4-1) 


where N_ is ak by seven deterministic matrix relating 6s to 8q. 
Ss is the state to be computed 
q is acolumn vector whose k components are the meas- 


urement data. 


The matrix N may be derived directly from the state transition 
equation written between the time for which the estimate is desired, 


t = ie and any previous time, t = t,. 
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[ 


rl 
_ =i 
is. - Tos -A, a, A.B &a dt (4-2) 


1 


where TS Se ee Nees For simplicity, the integral term in equaiiem 


(4-2) may be considered as a vector. 


i AG | 
sei) n a 
ae i yo Oa dt = Av ? (4-3) 
Ui Awa i 


and the transition matrix a asa partitioned matrix of partial dérivas 


Vvccwe Edict! Only the position vector ian t is considered, one obtains 


Cia oi Orne, fe 
a bv (4-4) 
Lin 
n 


where the matrix of partial derivatives consists of the first three rows 





oe 





or Ov om 
tl i) 


of Tn Hoquation (4-4) may be writlen for any nlimbemter Mime. 

t. GQ=n-i1,n-2....). An analogous equation may be written tes 
the propellant measurement and the results arranged in the tornigen 
equation (4-1). One might expect the propellant measurement to be 
less critical and less subject to error than the other measurements. 
Consequently it may be necessary to include only a few measurements 


of propellant, perhaps one or two. 


Let us consider now the problem of finding an unbiased filter, Fo: 
which minimizes the error in the state vector and satisfies the rela- 
tion 


68 = F.8q (4-5) 


If the uncertainty in the state vector is u and the error vector associat- 


ed with measurement is 


Mm 


, then 


6q=6qte = N6s+ Nu (4-6) 


}aQ zd || 


Since the true measurements and the actual state vector satisfy equa- 
tion (4-1), then 


= INU (4-7) 


The equation for the error ellipsoid associated with the measure- 


men, vector is given by 


eee 4 (4-8) 


where 


ee (4-9) 


is a k by k matrix and the brackets denote the mean value. Equations 
(4-6), (4-7), and (4-8) may be combined to yield | 


Mee 1 As aly 


vw N’ E“ Nu=(6q" - 6s" N')E* (6q - Nés) (4-10) 


Ij tne partial derivative of equation (4-10) with respect to the compon- 


ents of the staie vector are set equal to zero, one obtains: 


ie 
ou s 
2{—) nh’ & 1 Nu=-2N! BE”! (5g - Nés) = 0 
9 — =t = (4-11) 
Ss. 
i 
The solution to equation (4-11) is the estimate, 68. 
a -] . ae 
69-(N’ EN) NE! &q (4-12) 
Thus for the unbiased filter. equation (4-12) yields: 
2 =a 5 
a a ee Ne bee as) 
It is not surprising that the filter of equation (4-13) is the same filter 
a2 


derived by Potter and Stern using the method of maximum likelihood. 
4.4 Biased Estimate of Present State 


The proof that a biased estimate will result in a smaller ellipsoid 
oi error will be omiited since this topic is well covered in the litera- 


30,31,32 


vials In this section only the method of obtaining the biased 


filter from the unbiased filter will be presented. 


The optimum biased filter, Pop: is obtained by deleting the last 
seven columns from an associated unbiased filter which is computed 
using afictitious measurement of present state. Thus the measure- 
ment Crror VECtOr lo WSe [or computationwor the error covariance 


Matrix E is 


al 


€ = a (4-14) 


- bs 
The last seven rows of the N matrix must then be the seven by seven 
MGeTLIity Matrix, DNe reswlyiire sinters ecient en Fo computed from 


equation (4-135) but with the last seven columns deleted deiines tlie 


optimum biased filter, FoR 


45 Covariance Matrix of Neasnrement Brror 


The matrix Eis extremely important in the computation of ide 
State estimate. Since it involves the measurement errors trom seus 
discrete and continuous measurements, it is worthy of closer examina- 


aleval 


In section 4.3 a subvector of the measurement vector oq was 


written in the form 


Oe Oe (4-15a) 
From equation (4-3) it is apparent that 
n 


ow _ -] ~ 2 
eee | ep OOr aye J AB 6a at (4-15) 


where los the three by three identity, O-mis7a tinee byay ieee sume 


3 
iiave ix ane O 


3 
a three component null vector. Thus a measurement 


of position ee of two parts: the celestial fix, br, at time t = te, 
and an additional vector Ar, which contains the integrated engine vari- 
ations from t to t.. As a consequence, the error vector at time t. 
will also consist of two parts, 

t 


li 
: =f . 
é,= £4 +01, 0, O31] A; J ABe dt (4=1en 


1 


which for simplicity will be written as 


O2 


Oveewee |). Dic di (4-17) 


mmere € ee is the error in determining position and Co is he error all 


measuring the engine quantity. 


It appears justifiable to assume that measurement errors of dif- 
Beet types, i.e. position, propellant flow, acceleration, are indepen- 
dent and have zero means. Further, that the position measurement 
errors for two different times are independent. The assumption of a 
Zeomean ior measurement errors does not result in a loss of gener- 


ality, however it greatly simplifies the mathematics. 


With the above assumptions, the E matrix includes terms of the 


following form: 


1) On the diagonal 


i ls 
Tare | ie dt f €+ D- at> (4-18) 
See = | — a i. 
1 1 
2) Off the diagonal 
t t 
gl i 
< f De. at f ec! dD? at> (4-19) 
i ° ape Somes 
‘ J 


Using the theory for handling random processes”” Lae Inbar! 
terms are easily reduced to the form 
t i 
n n T | 
oie) (tyty) D> (ty) dt dt, (4-20) 
t: lie 
: J 
where @ (t,t,) is a diagonal matrix of autocorrelation functions. To 
proceed it is necessary to make some assumption about @. The most 
easily justified assumption is that the engine measurement errors 
each contains much higher frequency components. thandoes the matrix 
D, and each is uncorrelated except over short intervals. Thus, over 


the period of integration the measurement errors approach a ‘'white 
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noise’ distribution. With this assumption @ is reduced to a diagonal 
matrix of constants representing the measurement error variance 


multiplied by the Dirac delta function. 


o* 0 
G22 eae 8 (to -t.) (4-21) 
0 of eee 


(Measurement of propellant flow and exhaust velocity are used as ex- 


amples in (4-21) ) Designate, for simplicity 


6 = ), Sig te) (4-29) 
where y is a diagonal matrix of measurement variance. 


One integration of the terms in (4-18) and (4-19) may be performed 
to yield 
ba 7 
(Dee eas ct (4-23) 
Up, J 
where the interval k to n represents the shorter oi the intervals items 
or jton. The E matrix formed from these elemenis is a symmetpm 


positive definite matrix. 


The necessity for including engine measurements in the estimate 
of state for low-thrust vehicles is made evident by equation (4-24). 
Note that the measurement error vector, when position data alone is 
used, would contain integrals of the actual engine variations. That is 
t 


nh 


a J D, oa dt (4-24) 


1 

Assuming that the engine measurement errors are much smaller than 
the engine variations to be measured, equation (4-24) would result ina 
sizeable increase in the error ellipsoid. In addition, as the interval 

1 to n increases, the value of the celestial fix at time te is rapidly de- 

graded due to the ecifect of the increasing value oi the integralieriiam 


the HE mairix. 
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4.6 Estimate of Future State 


Ps snmeved previously, the future state vector of a continuous-thrust 
vehicle is not a deterministic function of the state history due to short 
period random variations and, perhaps, long term degradation of engine 
performance. If engine degradation, presumably in the form of loss of 
specific impulse, does not become a factor, the space navigator would 
be justified in assuming a zero mean for engine variations. With such 
an assumption a filter, biased to include a priori information of engine 
statistics and a priori information on state statistics, may be construct- 


ed which will produce an optimum prediction of future state. 


If long term variation in engine performance is evident, the naviga- 
tor is faced with the problem of extrapolating accumulated engine data 
ieerder tO make a prediction with acceptable confidence. Certainly, 
sophisticated analytic techniques exist for smoothing and extrapolating 
the measured engine data. However, in this section we shall only be 
concerned with using data, regardless of the manner in which they are 
processed, to predict a future state. For example, simple graphical 


extrapolation of plotted data will suffice 


Two prediction techniques will be presented, 1) a very simple ex- 
trapolation technique for short term prediction and 2) an optimum filter 
fecumique tor long term prediction. The first method assumes the 
present state has been determined by the methods of section 4 3 or 4.4 
ieuises the state transition equation and an estimate of engine perform- 


pee tO Predict State in the near future The estimate is 


tn+1 7 
A Z A aol A 
Se eis nr en | Oe oat (4-25) 
tn 
where the notation t = Eat will indicate the future. This method may 


be characterized as an extrapolation of filtered data and may be quickly 


computed. 


The second method uses the optimum filter of section 4.3 but in 


this case a representative subvector of the measurement vector is 


OO 


Sei seca | Dy Sa et < D, 64 at (4-26) 
1 n 


.th ae . 
the error vector associated with the i position measurement is 


ty tn+1 a 
: ot D,€, dt + J D.€, i (4-27) 
i 


n 


Ce. 
ak 


where 62 and €. represent the extrapolation of engine performance and 


ime error in thatecxtrapolanlionw: capectivel,. 


The optimum biased filter of section 4.4 may be used without al- 
teration; however, the covariance matrix of measurement error will be 


GE is hua Ulva iles Crore youre: 


This second method may be characterized as a filtering of extrage 


olated data and should be used for long range prediction. ‘The proof 


Dy otte reand Stern” shows this second method to be optimum and 


intuitively it appears to be of correct form. By comparing expressions 


for the covariance matrix of uncertainty 14 state at time toed tor thie 


tb 
two methods. 


U=<uul> (4-28) 


it is apparent that both methods result in the same error ellipsoid as 


tat approaches ty For both methods: 


El 
u=(N' BN) (4-29) 


t — 


lace 


thus the assertion that the simple technique will be quite adequate for 


Short term predictions appears justified. 
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Chie Pei vy 


TRAJECTORY DETERMINATION 


eet ouMmmary of Chapter V 


In this chapter the problem of finding a trajectory which satisfies 
the end condition, the propulsion restriction and which minimizes pro- 
pellant consumption is discussed and illustrated with the VSI trajectory. 
wie problem is formulated using the classical methods oi the calculus 
of variations. A new interpretation is presented for the results. which 
are shown to yield the same control program as the Pontryagin maxi- 
Mmieaeprinciple. The guidance theory of Chapter III is therefore sug- 
gesied as a trajectory computation scheme. Finally, qualitative as- 


pects of low-thrust optimal trajectories are discussed. 


mec ceneral Remares 

A vast amount of effort has been directed in recent years toward 
itesstidy Of Optimization problems. Such problems belong to the cal- 
culus of variations which owes its early development to such men as 
Lagrange, Euler, Hamilton and Gauss. The introduction of high speed 
computers has been the primary impetus in bringing renewed interest, 
after years of limited application, to variational Tec Ane le come Cia. 
ealy, much recent literature treats the characteristics of, the neces- 
Sary conditions for, and methods of computing solutions to optimization 
problems. Variational techniques are used almost exclusively as the 
primary mathematical tool. The contributions listed as references in 
this report constitute only a minor fraction of the published works. 
These efforts have resulted in a large body of theory now called opiti- 
mali control theory which has application to virtually all optimization 


problems dealing with dynamical systems. 


A fundamental precept of the theory is that along optimal trajec- 
tories, admissible first order variations in an unconstrained control 


Proctam Canney DROGHce atts order cil@el In tie cost tnciion, ~Hor 


of 


unconstrained controls, admissible control variations are those which 
do not cause a first order change in boundary conditions. This princi- 
ple is derived directly from basic equations and is part of the definition 
of optimality in the calculus of variations. A mathematical consequence 
of this principle is that certain arrays of coefficients, when evaluated 
along the optimal trajectory, will have a zero determinant. When work- 
ing with the practical problem of computing trajectories, or of guidance 
around the optimal trajectory, the inverse of a matrix defined by such 
arrays usually appears in the equations. Obviously the occurrence of 

a Singularity complicates the procedure of finding an optimal control 
program and its associated optimal trajectory. It has led to investiga- 
tion of second variations of the cost and the ea eCror ae to various 
schemes for fiidine tear -Opular controis!?, and to optimal controls 


which only approach the desired boundary conditions** 


The implication of these investigations is disturbing from a physi- 
cal viewpoint since they imply that an otherwise well behaved, smoothly 
operating system, in some sense becomes uncontrollable along an 
OpiUiIMalsreleuence (eaecror yal reality, singularities are more often 
mathematical than physical. A further implication is, that although a 
control can be found which approaches the optimal conirol to within a 
small increment, it is orders of magnitude more difficult to find the 


exacr Opilmal conirol, 


The physical world does not usually behave in such an unruly man- 
ner, thus the answer must be: 1) the mathematics have an interpreta- 
tion that has been overlooked or 2) the problems may be approached 


from a different viewpoint. 


Actually, both 1) and 2) have validity. To support this contention, 
the problem of generating an optimal reference trajectory for use wm 
testing the guidance theory of Chapter II] is considered as an example. 
To be sure, a conservative field such as the gravitational fieldis a 
well behaved space in which to work. Undoubtedly, problems which 
deal with dissipative forces or higher order nonlinearities may present 
difficulties not readily resolved by the methcd of this chapter. Hope- 
fully, however, it will provide an approach that can serve as a starting 


point for more difficult problems. 
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peo lhe Calculus of Variations Problem 


In Chapter II it is shown that minimizing the acceleration integral 
is equivalent to minimizing the propellant consumption for a low-thrust 
vehicle. Thus the mass rate equation is superfluous for VSI trajecto-- 
ries anditis necessary only to work with the equations of motion and 
mieeacceleration integral. For the VSI vehicle the problem can be for- 


mulated as: Given 


TNO are, (9-1) 
= (9-2) 
satisfy 
r(t.) =r. (33) 
V (te) a, (9-4) 
subject to 
ee (aie 2) 
Pade ae (5-6) 
: 
and minimize te 
I= if aM dt (5-7) 
2 0 


A scalar functional, F, will now be formed using the well known 


Lagrange multiplier technique. 
a? 

4 (5-8) 
2 


cL 


eee ane Z 
ee, 


(v + —z Fr - a) 
iP 

where A. and De are time varying Lagrange multiplier vectors (or 
Euler variables). The remaining variables have the same meanings 


as in previous chapters. 


mince the bracketed terms are zero, clearly 
Mi 
Nae emee ec (5-9) 
0 


og 


Integrating by parts and setting the first variation of the integral equal 


to zero, one obtains 


t 

6 f rat-o-] at eras by | ‘ 
0 

sy" Gere eae (5-10) 

= Sy Sie 


+ (a rf) 8a | at 


| @ 


Applying the fixed end conditions and setting the coefficients of the 
state and control variations equal to zero gives the Euler equations and 


the boundary conditions. 


\ va 
op + : 
| en Oo, L 
ee a * (5-1 
ES i Ay “G O3 
7 ‘\ 
ot = (5-12) 
{| —V 
\ } tf 


tr : a (5-13) 
| 


a=A (>=) 


Notice that the Euler variables, (5-11), satisfy the adjoint relationship 
to the variational equations of motion and that the boundary conditions 
are unknown. Equations (5-12) and (5-13) are written to emphasize 


the unknown boundary conditions. 


Clearly, if the correct initial value of the six component Euler 
vector were known, the entire system of state and Euler equations 
could be integrated simultaneously from t = 0 tot = te. The unique 
reference trajectory and its associated optimal control program would 
then be known. A procedure for finding the initial value is discussed 


in later sections. 
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Pevery annoying fact about the preceding approach to finding ex- 
tremals in the calculus of variations is the necessity of searching for a 
formulation of the problem which gives a meaningful result. Minimizing 
the cost function J is known to be equivalent to minimizing propellant 
consumption, thus one differential equation of state was eliminated, in 
Periicular, the differential equation describing the rate of change of the 
optimized quantity, mass. A simple answer resulted. It should not be 
Meeessary to search for an equivalent formulation for a problem if the 
Set O1 differential equations describing the system is complete, linearly 
independent and reasonably well behaved. ‘Two different Mayer formula- 
tions of the VSI problem help illustrate this argument. 

For the first problem, minimize propellant consumption instead of 
iienacceleration integral. Equations (5-1) through (5-6) hold, but the 
ivacs rate equation must be reintroduced, Thus it is required that 

m ear (5-15) 
2p 
and that ty 
m= 2 -m dt=m, - m, (5-16) 
be a minimum where m_ is propellant mass. For fixed initial mass 
Mie tirst variation of a equals the variation of - My. Horm in Ste 


mometional F 
ee 
ry 





FAL (e-v) tal (vt+-eer-ajtay (mt } =m 
r 2p (5-17) 
then integrating by parts and setting the first variation of the integral 


equal to zero as before, yields 
i 
f 


= tr 
6 f Fat=0- [AT or+aTov+a ~ 1) 6m 
0 —r = =v = m 


0 


(5-18) 








In this case the Euler variables for position and velocity again satisfy 
equations (5-11) through (5-13). In addition 





= = °? i 
(A a : (5-19) 
eal (5-20) 
i 
Aen 
ee = 0 (5-21) 
mM mm p 
Xr 
ie ——5- (5-22n 
eee 
ep 


Equation (5-22) appears to differ from (5-14) by a scalar function of 
time, since the control programic uilgue (9-27) mucuecaucesre 


(5-14). By manipulating (5-15) and (5-22) it is possible to obtain 


—e 


ee (5-23) 
es m 
which has the solution 
r a ee (5-24) 
m 
Applying the boundary condition (5-20) 
2 
ae 
m 
Substituting into (5-22) ; 
eee oN = 
a ae Nese (5-26) 
f 
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Equation (5-26) differs from (5-14) only by a constant which is easily 
absorbed in the unknown boundary conditions on Nee The Euler variable 


don appears to have been superfluous. 


Consider a slightly different formulation which is just as valid. Use 


thrust as the control then 


fe (9-9) 
f 
vi a (9-27) 
r 
° °2 
m=er-f— (5-28) 
2p 
te 
minimize m, = ecard (5-29) 
0 


Proceeding as before, the functional F and its first variation are 
f Z 





F=A° (r-v)tal (Vtyr-—) tA (nu $= =m 
ee ae i 2p 
(5-30) 
i te tp 
6 f Pat=o=| 5 {( ) dr 
0 0 0 
Vette i 
+( Jov+ (A -A ) dm 
- m —V we 
(5230) 
a 
7 
+(e - TF of | dt 
m oe - 


The blank brackets contain the same terms as (5-18) and are used to 


show the differences between the two formulations. In this case 





f 
ae 
a ee oD a 

ey 

r 
f = ——_- (5-33) 
Sem 

we 
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Using (9-33) to eliminate A’ from (9-32) and using the state equation 
(5-28), it is found that again 


A. =-2A_— (5-205 
solution (5-25) holds thus 


f=-, ma, (5-34) 
ce 


which is the same as (5-26) since a = f/m. 


The preceding VSI examples result in the same control and in none 
of them is it necessary to use the variable Aen Betore interpret, 
what this means one additional example is presented which involves 
placing a constraint on the thrust. Assume that 


2 2 


mS (5-35) 


where Be is the thrust limit. Following Breakwell? ', the constraint 
may be handled by adding to the functional (5-30) the term 
! 
y (f° - f<) = 0 (5-36) 


J 
where y is selected to satisfy (5-36). Thus for this case 


F={(5-30)}+ ye f“) (5-37) 
ve te — te 
6 f Fa-o-| | a [ Ore ) Ov 
0 0 0 
T T 
| f . f 
{p= —V eae ! T 


(5-38) 
where the blank brackets are the same as (5-31) and not needed in the 


ee Culler, 


Consider the coefficients of 6m and 6f in (5-38) 


Miner = 0 (5-39) 


al 
m V 


Ay 
a : (5-40) 


| ) 
r =(y'+— ae (5-41) 
rr 


When f # f y' must be zero to satisfy (5-36). Therefore the 
solutions for Ne and f in such unconstrained regions can differ from 
(5-25) and (5-34) at most by a constant. When f = i the only admissi- 
ble solution for (5-40) is 








f 2X 
= eee (5-42) 

~ fal 

— V 
Therefore 
Br: Ir, | 
€ +2) ee (5-43) 
p ie 


With this result it is possible to rewrite (5-40) with a new variable y 


een that 


f=mya (5-44) 


where 


(5-45) 
mld | 
—V 
The Euler variable Ne and its differential equation have again been 
eliminated from the problem. 


In the preceding examples the manipulations required to remove 
An are not particularly difficult but by no means are they obvious. ‘The 
meaning of a superfluous Euler variable immediately comes into ques- 


tion and it seems reasonable to inquire if there is general significance 
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in this phenomenon or if the preceding examples are only isolated 
special cases. Conceptually there is nothing about the variable A 
which is unique except for its association with the cost variable m in 

a Mayer form problemy In reference 17 BGreakwell inmierprers ome 
Euler variables as the sensitivity of the absolute minimum of the cost 
function (i.e. the unconstrained minimum) to changes in the state vari- 


ables. Therefore ior ine proplemmvarsiis tests 


dm am. 
At (t) = p = . f 
Os, Os, 








(5-46) 


where m”* denotes the minimum value of m_ obtainable from the initial 


state when no control constraints are used; similarly for - Me - 


In a typical two point boundary value problem with physical rates 
of the form s = g (s,a), the usual objective is to find a control solution 
which satisfies the end point. If any of the state variables have free 
end conditions the solution is not unique. Optimization criterion are 
then used to assure uniqueness. From Breakwell's interpretation, 
when the Euler variables are formed iric a vector they describe the 
direction in state space of maximum sensitivity of the cost, i.e. the 
gradient. However for state variables involved in the cost function 
such information is available directly from the governing differential 
equations. In answering the question: ''How do state variations affect 
the cost?"', the Euler variables provide a coupling between state vari- 
ables and the cost function which for some variables is not evident 
from the system equations. But for state variables involved in the 
cost this coupling is essentially a priori information and is in some 


sense superfluous. 


The classical approach in the calculus of variations assigns an 
Euler variable for each state variable. If a variable is superfluous it 
should drop out of the formulation. But finding a way to assure that it 
does may be exceedingly difficult, even when it is recognized that the 
variable is superfluous. As is subsequently shown, the Euler variables 
correspond to the costate variables in Pontryagin formulations. If the 


variable is superfluous in the calculus of variations it is also superfluous 
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using the Hamiltonian approach. This fact was tacitly acknowledged in 
Chapter III when the state variation, om,,was dropped from considera- 
tion and the last row of the adjoint set was deleted. The reason given 

in Chapter III for deleting one row of adjoint functions was different 


from here, but the result is the same. 


The general significance of a superfluous Euler variable is indi- 
cated by comparing the classical approach to the newer state space 
approach using the Pontryagin maximum principle The classical 
formulation suffers from several difficulties, two of which are: 1) 
The number of variables and equations are usually quite large and. 
perhaps due to tradition, are often treated individually as scalars As 
a result algebraic detail often obscures important ideas in classical 
formulations. The more compact vector and matrix notation are only 
now becoming widespread in the literature. 2) The theory offers no 
Suggestion of how to solve for the correct boundary conditions on the 


Euler variables. 


However, an important attribute of the more detailed notation is 
indicated by the preceding examples’ By sufficient manipulation an 
unneeded variable was eliminated. It is more than coincidence that 
the variable is precisely the one that creates problems in the state 


Space approach. 


The state space approach using vector and matrix notation as in 
Chapter III is popular because of its compactness. Further, when 
coupled with the Pontryagin maximum principle it may be used to 
solve optimization problems including cases of linear control (the CSI 
problem) which can not be completely solved in the calculus of varia- 
tions. The approach also provides a method of determining the bound- 
ary values for the costate variables. Finally, the optimality criterion 


is more useful. This last statement is discussed fully in Appendix D. 


In the state space approach the system differential equations are 
written in vector form, an appropriate cost function is chosen and then 
methods are sought to find the desired solution. From the references 
previously mentioned it is apparent that the solutions often encounter 


a singular matrix. In the problem of this thesis the singular matrix 


67 





1s denoted as M, and is the seven by seven matrix corresponding to 

M in Chapter Ii]. These matrices are discussed in detail in Appendix C, 
It is suggested that such singular matrices may be the result of attempt- 
ing to optimize a function of state and at the same time retain a costate 
variable (or Euler variable) which is superfluous. In the problem of 
this thesis such was the case. li is further suggested that difficulties 
sometimes encountered in numerically finding terminal values for the ; 
Euler variables may also be due to the presence of a superfluous Euler 


variable which was not recognized as such. 


No attempt will be made in this report to derive the general rule 
which will show when an Euler variable is superfluous and thus may be 
eliminated a priori. It may be the case that those matrix elemenm 
which cause nonphysical singularities to occur in state space formula- 


tions can always be removed without changing the problem. 


In the remainder of this chapter the relationship of the calculus of 
variations solution to the state space solution is shown for the low - 


thrust problem only. The method of deleting Nes is illustrated. 
9.4 Removal of a Superfluous Euler Variable 


Although an Euler variable may be superfluous, a general method 
of removing one is not obvious from the derivation in section 5.3. Con- | 
sider, however, the following derivation of the problem in section 9.3 
which combines the compact notation with the classical approach. It 


may be indicative of a general method. 


Assume the system is described by physical rates which have the 


LOrim 
s = g(s, a) (5-47) 


where s is the state vector, a the control vector and g is a vector funigs 


tion Of State and control. 


Assume that the cost variable is some scalar function of the physi- 


cal rates or some nonlinear function of the control. 


Se ana) (5-48) 
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or 
Sia) (5-49) 


where S is the cost integrand 


Further, allow the physical rates to be subjected to constraints. 
Define the constrained variables as a vector function of the physical 


mates or of the control. That is 


ie ees eo (9-90) 


A general functional, equivalent to those used in section 5.3 is 


F=Statgty'ly (5-51) 


—_ 


where A is the vector of Euler variables and where Me is determined 
in such a fashion that the constraint is satisfied It is a vector equiva- 


tem, of the v in section 5, 3 


Application of the usual variational technique produces the expanded 


set of state equations, Kuler equations and control equations. 


Ss = g (s, a) (5-47) 
. aS .T a ‘Tay \t 
hei Ds SB +y —X | (5-52) 
- Os = Os as 
as T 28 rE ay 
o- (2 +a °B +y i (5-53) 
_ da 7 da da 


For the unrestricted case such that ca O and using the familiar 


Melimitionol the mairices A and B, that is 


OB A (5-54) 
os 
oe (eeaen 
da 
mien from (5=52) one obtains 
ee (5-56) 


De 


A solution to the set (5-53) to within an unknown scale factor, is 
a=B' (5-57) 


That (5-57) is a solution may be verified by expansion. The variational 


equations of state and their adjoint set are 


bs = A 6s + Biéa (5-58) 


A=-AA (5-59) 


Since the Euler equations (5-56) satisfy the adjoint relation as does 


(9-99) a solution for (5-56) is obtained by choosing 
Nail (5-60) 
Then 


: ak - 
Ayr AD dy (5-61) 


Dropping the subscript t from (5-61) and substituting into (5-57) yields 


| © 


=B AX, (5-62) 


By partitioning /\ and r¢ it is apparent that deleting the seventh row 
ola US iy and deleting the unknown boundary value Ame eliminates 
the superfluous Euler variable from consideration. That is 


rv 
—r 


a-BT/A*™! Auld Ay (5-63) 


In order to find the optimal trajectory it is only necessary to fima 
the final value of the six Euler variables. (Or, since the adjoint and 
fundamental solutions, discussed in Chapter III, allow transformation 
at will between one terminus and the other: the initial values of the 


Fuler variables may be used. ) 


By letting 


ee ae (5-64) 
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and deleting eae and my from (5-63), the expression yields 


a=B! A*? py (5-65) 


which is identical with the expression for the optimal control program 
for guidance derived in Chapter IJ]. Thus beth the Pontryagin principle 
and the Euler variables give the same result when he is deleted. Ifis 
showr in Appendix C that along the cptimal trajectcry A is a null vec- 
ror of Bt A thus equation (5-62) displays the wel! known and very 
troublesome singularity when evaluated alcng the cptimal trajectory. 
This is the consequence of optimality which has instigated the search 
tor new methods,including this cne. It is further shown that V1: net a 


: a ace. | 
null vector cof BO A thus (5 65) can be used. 


For the restricted VSI case the preceding arguments hold with y 
Savering as in section 5.3. The restricted VS! case is completely 
determinate in the calculus cf variations using Breakwel) s approach 
and deleting a The CSI case may alsa be trea’ed 1n an anaicgeus 
fashicn except that the ccast phase can nc’ de uniquely determined. 
This 1s due to the so-calied degeneracy of the calculus of variaticns 


for linear optimums 
9.9 Solution by Direct Integration 


There are two general appreaches fcr finding an optimal trajectory. 
Descriptions of the first methcd are usually prefaced by: "If anv non 
optimal solutien can be found which satisfies the boundary conditicns. 
then the solution can be mcved in the direction of the optimal, etc." 
This approach is usually called 2 gradient methed. It 1s rot used in 
this report because nonoptima! sclutions carnot satisfy the Euler 
Equations, thus it would be difficult te relate any nonoptimal ccnirci 


to the form of equation (5-65) 


The second general approach is often described as ‘'sclving the 
wrong problem in an optrma! manner '' This 1s the approach used in 
Chapter Vi. An initial guess 1s made for the initial Euler variables 
V. then the state equations, adjoint set and Euler variables are integrat- 
ed to the final time. The resulting final values are compared with the 


desired final values and a new estimate of Vis cempuited. The precess 
t p 


cA 


is repeated until the boundary values are satisfied. The fact that the 
gravitational field is conservative and well behaved aids materially in 


improving the Speed ol Convercence. 
5.6 Properties of Optimal Trajectories 


It is possible to gain considerable insight into the properties of 
propellant-optimal trajectories without resorting to machine computa- 
tion. Such insight aids materially in working with the mathematics of 


the transfer process and often leads to new methods of approach. 


For this purpose, consider the problem of transferring betweentwo 
planets with minimum propellant expenditure as one of changing energy 


and angular momentum in the most efficient manner. 


From the energy integral of orbital mechanics one may write the 


energy per unit mass of the vehicle as 


2 
Pee ae aye (5-67) 
y ig 
and the enerey lime derivative as 
St +Fo x (5-68) 
ic 


Equation (5-68) may be written in the vector notation and reduced with 
help of the equations of motion to 


T 


era (5-69) 


I< 


Equation (5-69) represents the rate of energy change imparted to Gae 
vehicle. The rate of energy expenditure by the propulsion system is 
the power 

aes ae 

ae wef (5-70) 


foe + be 2 





where m is considered as the mass flow rate per unit mass and c has 


the opposite sense of a. 
similarly a functional expression for angular momentum is 


h=rXv (5-71) 


eg 





herXvtrxvy (5-72) 


where the notation signifies the vector product of r andv. Equation 


(9-72) 1S reduced with help of the equations of motions to 
h=rXa (5-73) 


Using equations (5-69) and (5-73) the problem is formulated as: 


t 
f 
Minimize ST Sea ch (5-74) 
0 
subject to the constraints 
tp 7 
Wei ch (5-75) 
0 
My 
Ah= f{ rXadi (5-76) 
ae 


plus velocity and positicn constraints and where Ae and Ah are known 


quantities. 


The above formulation is not computationally desirable’ however 


it provides a basis for some significant qualitative arguments. 


Jt is apparent from equation (5-69) that with a fixed in magnitude 
the vehicle changes energy most rapidly, thus most efficiently, when 
vis large and v and a are colinear. From the energy equation (5-67) 
one observes that in a central force field the velocity magnitude goes 
inversely with the square rect of the radius From equation (5-73) it 
is apparent that h changes most rapidly when r is large and orthogonal 


to a 


Fram the preceding cbservations one may deduce that the trajec- 
tory which minimizes the propellant expenditure will tend to keep the 
exhaust velocity vector aligned with the large velocity vector when the 
vehicle is deep ina gravitational field so that energy is changed most 
efficiently. It will also tend to rotate the angular momentum vector 
when the vehicle is far out 1n the gravitational field sa that the vector 


r has large magnitude. Where these requiremerts are incompatable 


ho 


with boundary conditions, the optimum solution should resolve the 


conflict in favor of propellant conser, anon 


The trajectories which have been generated numerically in this 
thesis all appear to satisfy these precepts in so far as boundary con- 
ditions permit. The increase in specific impulse which is character- 
istic of variable thrust rockets as they traverse the center portionug@s: 
the heliocentric phase is due to the acceleration vector becoming ofuiaiag 
onal to the velocity vector. Likewise the coast phase for a thrust- 
limited rocket occurs when the acceleration vector rotates throwals 


the orthogonal orientation. 


With the qualitative arguments of this section, it is possible to 
sketch a reasonable approximation to a low-thrust optimum trajectory 


without resorting to machine computation. 
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COMPUTATIONAL PROCEDURE FOR TRAJECTORIES 


6 1 Summary of Chapter VI 


Details of a procedure for using the guidance equaticn as the basis 
meer eajectory computations are set forth in this chapter. Specifically. 
miemcaiierential equations to be integrated and the procedure for cor- 


recting initial conditions on the Euler equations are presented 
pee General Remarks 


In Chapter V the fundamental problem of selecting initial (or final) 
mem@eriions on the Euler variables is intraduced Since the acceleration 
program is an explicit function of the Euler variables it is necessary to 
meke a tirst estimate of the Euler initial conditions in order to start 
the iteration ‘This estimate will usually be a gross error under the 
est Oi Circumstances. thus it 1s desirable to choose a very simple 
estimate such as the null vector or a unit vector. Because of this, the 
mpeoccdure we seek must be strongly convergent and independent of the 


error in the first iteration attempt. 


iietiay erative procedure tie speed of Conver sence is directly 

dependent upon the validity of assuming that certain quantities do not 
change from one iteration to the next. In the preceding chapters it was 
explicitly assumed that the miss vector at the terminal point was de- 
pendent only upon a change in the initial Euler vector and that A and 

B were invariant for neighboring trajectories. For small perturbations 
around a reference trajectory such an assumption is valid. However 
for the large perturbations expected in trajectory computation i \aena 

B are not invariant. The change of these quantities from one iteration 
feoernie mext and the change of other second order quantities will deter- 


mone tie Speed Of Convergence Of the procedure. 


eS 


ion 3 The Correction Procedure 


From Chapter III one may write the optimal acceleration program 
as a scalar times a matrix function of the adjoint set and the initiali- 


ZAnloOW Vector. 
ae 7 ve Bt Ax*xt ve (6-1) 


where y is the scalar switching function derived from the limiting con- 
dition on acceleration, Bis the seven by three matrix of partial deriva- 
tives, 9g/8a, A_* is the six by seven reduced adjoint set and v is the 


six by one initialization vector. 


In order to simplify the notation in certain equations to follow, 
Euler variables will be used interchangeably with the matrix notation. 
Recall from Chapter V that 


-=B  A*” v (6-2) 


and 
rn, =- 7. (6-3) 


= 


The objective of the iterative search for the optimal accelerates 
Prog tai WS to. itidetie aycelor vo which causes the six component miss 


vector, €, to approach the null vector to some desired accuracy. iaee 


ee, ay (6-4) 


—n n-l 
and we desire that 


ao 7 
Ue (Oy 


for n as small as possible, where n is the iteration number. 


The miss vector, €, may be written as the difference in final state 
which results when a nonoptimal acceleration program, a, is used in 


place of a. Define 


i 
f 
n= {| A*Badt = My (6-6) 
Bra Me 
-&= f{ A*B (a°-a) dt (6-8) 
0 
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then using (6-1) and (6-8) 
Ve 
-§:fA*BB A*! aiyy)jat (6-9) 
0 


Equation (6-9) implicitly assumes that JA. and B are invariant functions 
of time and thus constant from one iteration to the next. The assump- 


tion permits adequate convergence and is used throughout. 


itimorder to provide for the discontinuities which occur in certain 
variables at the beginning and end of restricted thrust regions, ii is 


menvenient to write equaticn (6-9) as 


j-1 ; tit] 
oS [A BBL A«t yu dt Av (6-10) 


where j is the number of distinct phases of the trajectory Expanding 


i100) according to Leibniz’s rule yields 


j-1 tit] : 
ee), f(A® BB At * y+ A* BB A*®* y%)at 
1=0 tj 


ee 


t. 
“Itl ot eel) 
aa al aet yYv “ Av 





if : 


i OV 


(k= 4, 141) 


> 


since y is the discontinuous switching parameter, the time at whicha 


discontinuity occurs, t, and its derivative with respect to v may be eval- 


k 
uated from the derivation in Chapter 1V. Rewriting equation (3-76) we 


obtain 
oq 
Ot Ov 
cs = glee (6-12) 
OV dq 
i ot 
By comparing (3-52) and (6-2) observe that 
= | Ge 
a= |r. | (6-13) 


ror 


Using (6-2) and (6-3) and taking the derivative, yields 


ot. ho Bi Ax? 
= = - =. (6-14) 
r 

— —v -r 


The term 8y/8v may be evaluated from the definitions in Chapter III, 


which are written here in more convenient form. 














a 
eee ele : | (VSI) (6-15) 
ee 
Q 
Ol ya mo l = il 
la | 
y = (CSI) 
2 
2 te ip 
SP sel en ee 
i ane (6-16) 


Therefore by differentiating and applying equation (6-2) 
a ho Bone 


3 ify # 9, 1 
In| 


= 0 Tie Ut 
(6-17) 


From equation (6-6) let us make the following definitions for 7 and M. 


eee BrAxTy (6-19) 


n= Mv (6-20) 


Using (6-11) and (6-14) through (6-20) one obtains 


oe too oom Vout off = 
YHA 7) 
a. i M - aoe 2 dt Av - T Av 
0 (2) vyIA- a. | 

a [eon 

(6-21) 
where we let | 0 oe 1, o| 

> (6-22) 
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The final term in equation (6-21) is evaluated at points such that y is 


not zero, i.e. € before cut-off ande after cut-on. 


Define 
t, oop ot ‘off 
= a ep dt - TF (Gi- 25) 
0 ( 2 | YIA 
—— —v -r 
m ‘ 
on 
then 
_ 4 
BU ee rc (6-24) 
and from (6-4), (6-5) and (6-6) 
vo=>Min Me te 16-25) 
minus the correction procedure to be emplayed is 
eae = me + S (6-26) 


= et 
6.4 Quantities to be Computed 


In order to obtain the quantities M”™ and € of equation (6-26) it is 
necessary to integrate the state variables, the adjoint set, the vector 


7] and the matrix M*. I[n additicn 1+ 1s convenient ‘a integrate the Euler 


variables rather than compute them from other quantities. 


From the discussion in Chapter V it may be roted that the adjoint 
equations are to be integrated backward in *1me fircm ie =]. However, 
if 1s easily shown that integrating backwards in time is an unnecessary 
complication. In Chapier IJ] the variaticnal equaticns of state and their 


adjoint set were solved together tc produce 


t 


f 
A;ts,2A,o,+ f AB tad 6-27 


t 


and BA was arbitrarily chosen equal to the identity in the guidance 
problem. For computing the entire trajectery from t=0 tot = t, one 


may write an equivalent equation 


+ 
oi 
se Ke 3c i: Ps T wT 
See Oe ee BOR May pec 
0 (6-28) 


fee) 


where A,” has both the last row and last column of A deleted in order 
to be compatable with € and where OS = Osince the problem starts 
from a given state. Since the initial condition on A. is still arbitrary 
it is convenient to choose N., = ] so that backward integration is not 


necessary. ‘Thus in place of (6-24) use 
-AS* E=M™ Ay (6-29) 
N=! (6-30) 


An interpretation of this procedure is as follows. An arbitrary accel- 
eration program is selected (1.e. estimate vy); then the state andwags 
joint equations are integrated forward resulting in a miss in position 
and velocity, - €. Since the acceleration program satisfies the condi- 


tions for an optimum to whatever terminal point it reached, namely 


aye 


Thus we are assured that when &-»O, m,.—»m also. It is un- 
~. — “ft ta (etetaNel oe ) 


important, therefore, that the ''miss" in m, is unknown, since if the 


is 
( oe ‘| , the propellant required to reach that point is minimum. 


procedure converges on the physical boundary conditions it does so via 


a propellant-optimal trajectory. 


ester im a € merely reflects the physical miss back to the 
initial state and compensates for the procedure of integrating the ad- 
joint set forward in time instead of backward in time. In the machine 
procedure it is convenient to use equation (6-31) in lieu of (6-26). 


~ ~ 1 de oie 3 
De ee gee a (6-31) 


The full set of equations that must be integrated are listed with the 


appropriate initial conditions. 


S = g (S, a) (6-32) 
SO) ge (6-33) 
A- -/A A (6-34) 
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ee 1 (6-35) 


O 
eee ry (6-36) 
A (0) =v (6-37) 
yo 
; a she sie Y tT) 
a ee (6-38) 
(<2) 
iy) 
M* (0) = O (6-39) 
n=Mv (6-40) 
where M is the first term of M” 
M (0) =O (6-41) 
a ane (6-42) 


An information flow chart and the FORTRAN program used to mecha- 


mize the precedingequations are included in Appendix G. 


The procedure worked satisfactorily. converging to very small 
values for § with surprisingly few iterations. The initial estimates 


for v that were used are 


De®) (Woltrarecrories} (6-43) 


(CSI trajectories) (6-44) 


Mo IO 


= 
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CHAPTER. vil 


A NUMERICAL EXAMPLE 


7.1 Summary of Chapter VII 





In this chapter a sample heliocentric transfer from Earth to Mars 
is described. The selected mission was programmed for machine com- 
putation to test the theories developed in the thesis. Numerical results 


from the test are presented and discussed. 
7.2 General Remarks 


The linear guidance theory of Chapter I]I assumes that an optimal 


reference trajectory is known. In order to test the guidance theory, 


. 


therefore, it is necessary to compute an optimal reference trajeciomes 

In Chapter V the discussion of a superfluous Euler equation shows that 

the solution of the linear guidance equation satisfies the necessary 

conditions for an optimum. Thus in theory, iterative solution of the 

guidance equation will generate the desired optimal reference. If the 

technique converges the linear guidance theory is validated. If the 

procedure converges rapidly then the theory also provides a simple, 


fast and effective way of generating low-thrust trajectories. 


The procedure was found to converge rapidly for both CSI and VsI 


vehicles. 
7.3 The Mission | 


The mission selected was a 150-day heliocentric transfer from | 
Earth to Mars. Flight time of 150 days was chosen in order that 
numerical results could be compared with a 162-day coplanar trajectory 
penerated iby Friedlander‘. The period chosen corresponds to the 
favorable opposition of Mars during the summer of 197 ot, The depame 
ture date, from a position on Earth's sphere of influence, is J. D: 244— 
1090.5. Arrival at Mars sphere of influence is 150 days later. The 
vehicle is assumed to have Earth's orbital velocity at departure and to 


match the Martian velocity at arrival. These velocity conditions are 
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mameopilinal for the entire mission, i.e. escape, transfer, capture, 
Mowever the author wanted to test out-of-plane components of the state 
vector. The use of planetary velocities for initial and final velocity 


Bermitied this. 


Both CSI and VSI modes of control were tested. The assumption 
was made that a space vehicle of given size and power could be con- 
mowed under either mode of operation, CsSl or Vsl. Thus two complete 
sets of data were generated which differ only in mode of control. In 
addition the adjoint sets for an optimal thrust program and an optimal 
feeeleration program were generated for each control mode. The pur- 
pose of having both types of adjoint functions is to permit all engine 
Paomlalies Ooi interest to be studied. The subtle difference between the 


two sets of adjoint functions 1s discussed in Appendix C. 
7.4 Computational Coordinates 


The coordinate system used for the mission analysis is a simple 
but effective system which is defined by the transfer plane. The trans- 
fer plane is the plane which passes through the sun and contains the 
weoited departure and arrival points. The x axis passes through the 
launch point, the z axis is normal to the plane in the northerly direction, 
wieetie Y axis completes the right hand triad, A procedure for trans- 
forming ephemeris data into computational coordinates is derived in 


mppendix F, Figure 7-a illustrates the geomeiry of the transfer. 
fo Engine Selection 


Engine sizing for the sample mission was computed on the basis of 


the mass distribution for maximum payload derived in Chapter Il That 


2 ee af (7-1) 
q Mo a | 


A trial trajectory for a mass independent VSI rocket was computed to 


is 





obtain a first approximation for J. The value for the one-way transfer 
was scaled up by a factor of 4 to provide a more realistic value for a 
round trip mission; also to insure that mass would never be reduced to 


zero during computer tests. A value for a of 10 kg/kw was selected as 
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being reasonable, though perhaps unattainable in less than a decade. 


The resulting value for exhaust power is 


oy TORO 24 ee (7-2) 
Mo kg 
This value for power was used for both modes of control, (VSI) and 


(CSI). 


The VSI trajectory was generated and the resulting value for initial 
feeelcration, 1.57 X i = go was used in the field-free space equa- 
tions of Appendix A to obtain an optimum initial acceleration (thrust 
level) for the CSI vehicle. A value of 1.2 X 10°" g was used. This 
fee Of Optimum initial acceleration, based on field-free space, is in 
surprisingly good agreement with the more rigorous computations of 


Melbourne and Sauer. 
7.6 Numerical Results 


The computer output data for the final iteration of each mode of 
control are reproduced as Appendix H. In addition, plots of several 
interesting output quantities are presented as Figures H-a through H-r. 
The data samples and the plots are strictly valid only for the particular 
meoes they represent, however they are indicative of the order of mag- 
nitudes applicable to many one-way missions in the solar system. The 
adjoint functions, for example, (Figures H-d through H-r) are in good 
agreement with the values obtained by Friedlander for a coplanar trans- 


fer. 


It is interesting to observe that the acceleration magnitude for a 
Wel vehicle, Figure H-a, is approximately a linear function of time and 
almost symmetric with respect to the midpoint. Comparison with the 
field-free space acceleration program, Figure 2-c, whichis linear and 
symmetric, confirms that to first order analysis of VSI trajectories in 


field-free space may be applied to the gravitational field. 


This result only confirms the work of other investigators and was 
aaticipated,. Ine analyvcis tor Col velicles imireld-tree space is not 
as straightforward and requires a large amount of work. The field- 


free space derivations in Appendix A result in a value for the optimum 
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CSI thrust level which is in good agreement with the computer results 
of Melbourne and Saver. The value obtained for the final mass ratio 
is larger than the value which results from computer studies in the 
gravitational field. Referring to Figure 2-d, observe that the point 
representing the CSI transfer falls slightly below the curve represent-~ 
ing field-free space prediction. This is satisfying in that the penalty 
for CSI control is less than expected. However, before concluding that 
field-free space is always a good predictor for CSI transfers, it will 
be necessary to test a large number of additional points to determine 
if agreement is sufficiently good to warrant the computational effort. 
The limitation of time has prevented the author from undertaking such 


a study. 


In Figures H-b and H-c the transfer plane components of position 
and acceleration are ploited. The out-of-plane components are quite 
small and are not considered further. A surprising result is that the 
maximum difference between the physical paths of the VSI and CSI 
trajectories is less than 0.01 A u. The large difference in the form 
of the acceleration programs, Figure H~a, leads one to anticipate 


rather large trajectory differences. This is found not to be the case. 


A very important characteristic of continuous-thrust rockets is 
shown in Figure H-p. This is the sensitivity of position error to mass 
change when the thrust program is specified. The y component of 
position at arrival time is subject to an error of 0.81 x oe km for 
1% variation in launch weight. For a one-ton vehicle this is approxi- 
mately equivalent io a 25,000 mile terminal error for each pound of 
launch weight variation. The large sensitivity to mass changes 
emphasizes the requirement for accelerometers with sensitiviltiesuom 
Ene woOrder OV 10 g. such instruments will be extremely important in 


the navigation of low-thrust spacecraft. 
(27% Convergence of iiescomputatouemnouLine 


The results discussed in the preceding section, although interest- 
ing, are only by-products of the experiment devised to test linear 
guidance theory as a computational method for trajectories. Conver- 


gence of the routine is crucial to the thesis since the linear theory 
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purports to null the miss vector for small variations in state. Appli- 
@amonm Of the guidance theory to the launch state, i.e. trajectory com- 
putation, provides the most severe test of guidance. The results of 
frevest are plotted in Figure /-b. ‘The criterion chosen as indicative 
of convergence is the RMS value of miss vector length. In order to be 
dimensionally compatable, the velocity components are multiplied by 
meat time, 150 days. Therefore the ordinate in Figure 7-bis the RMS 
So 


a: The abscissa is the iteration number. 
f 2 


value in A.u., of 
V 


eeveral cases were tested; the three plotted are indicative of all tests. 


dhe first one or two points in each case represent the result of 
estimating initial values for the Euler variables, v, and were widely 


Peanered, These points are not plotted. 


mie unrestricted VSl trajectory converges very rapidly; reducing 
|e | by approximately two orders of magnitude each iteration. Con- 
vergence is very smooth. The CSI cases converge less rapidly and 
less smoothly, but satisfactorily. In general, as the thrust limit, ao 
fercduced, convergence is slower. If, for example, ao is reduced € 
below the value which permits no coast period, the procedure will not 
converge. This physically corresponds to a vehicle "in extremis” 
such that insufficient thrust is available to complete the mission ‘The 
result also indicates that reserve power must be available for guidance. 
A. CSI vehicle which is beyond the coast phase cannot correct for state 
error in all six components of € unless reserve power 1s available. If 
reserve power is not available, guidance must be based on a formula- 


tion such as Pfeiffer's?*, which gives the minimum miss. 


It should be reported that the CSI miss vector magnitude entered 
a small limit cycle around the target prior to addition of the term 
ot / av. This term is derived In section 6.3. Prediction of the switch- 
ing time by including the above derivative is necessary to obtain the 


eesired accuracy, 


Ipimay be desigapile tonincinde an cdditional termi in tie preceamec 


to smooth convergence of the CSI routine. The term 9m/ dv, evaluated 
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at the switch points, will probably be sufficient to do this. Computer 
tests were completed before the significance of this term was recog- 


mezed and time limitation precluded retesting of the procedure. 


The accuracy obtained from the tests is significant. The routine 
was stopped after a predetermined number of iterations to conserve 
femipurer time, however for Vol trajectories, the largest error in a 
component of position was of the order of one mile. The largest error 
in velocity was of the order of 20 feet per hour. These values are given 
in the computer results in Appendix H. The computer was stopped 
weiore attaining this accuracy for CSI trajectories, however the con- 
vergence plots in Figure 7-b indicate such accuracy is attainable with 


a sufficient number of iterations. 


OW 


CHAPTER Vii 


SUMMARY AND CONCLUSIONS 


6 1 Summary 


In this report the obiective has been to justify the argument: ''There 
exists a linear method which produces a4 prcpel'ant-optimal control pro- 
gram in anoniterative form for guidance cf low-thrust space vehicles 
and which provides a simple, rapidly converging iterative technique for 


computing propellant-optimal trajectories, '' 


In pursuit of this objective 14 1s necessary, in Chapter Il, to develkas 
the parameters which characterize low-thrust vehicles. Due to ‘state 
of the art'' restrictions on variable-specific~impulse (VSI) machines, 
the parameters fcr constant-specific-impuise (CS1) vehicles are devel= 


oped an additicn 4o the 1deaiszed mode of ceniva. 


As an adjunct of Chapter II, the equations for field-free space appli- 
eaicn of CS! contre! are derived in Appendix A. This derivation is 
used to test the validity of approximating CSI transfers in the solar 
system by field-free space analysis. as has previously been done for 


constant-power vehicles. 


The parametric dérivaticns provide a starting point for investiga- 


ting propeliant-cptimal guidarce of low-thrust vehicles. For this study 


t i 


the vehicle 1s characterized by its ‘state The state is represented by 
a vec’or consisting of threé components af position, three components 
of velocity and the vehicle mass. The differential equaticn of state 1s 
linearized by considering variations of the state vector with respect to 
an optimal reference trajectory. Auxiliary functions, called adjoint 
variables, are used to solve the linearized differential equation. The 
solutior to the expanded se’. of equations is called the state transition 
equation or fundamental guidance equation and serves as the basis of 


the guidance theory. 
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The linear guidance theory is developed in Chapter III and applied 
to both modes of vehicle control. The solution for an unrestricted VslI 
vehicle is derived using the calculus of variations and then using Pontry- 
agin's maximum principle to prove that both methods lead to the same 
feoull. The Csi problem is then solved using the Hamiltonian. The 
Crucial step in either approach is deletion of the adjoint functions asso- 
ciated with the mass rate equation. Deletion of this set of functions to 
form a reduced adjoint set permits a solution to be obtained directly. 
Otherwise, the system is indeterminate due to a singular matrix. Proof 
of the singularity and justification for deleting the one set of adjoint 


Pmictions are treated in Appendix C and Chapter V respectively. 


As an adjunct of Chapter III, explanation of the adjoint relationship 
and derivation of Pontryagin's principle are presented in Appendices B 
and D respectively. Also in Appendix Dis a discussion of optimality 
criteria as derived from the Pontryagin principle and from the calculus 
of variations. IJmportant properties of the state transition equation are 


presented and discussed in Appendix C. 


ie problem of estimating ine yehicle state is studied in Chapter TV 
as a problem in navigation. The method of redundant measurements 
eeucied by Battin, Stern and Potter is extended to include continuous 
measurement of low-thrust engine performance. Based on the concept 
of filtering redundant data with a biased filter, two methods are derived 
Moteepredicting future state. One is a simple method for short term 
Prediction; the other is more complex but also more accurate for long 
Bese prediction. The effect of omiiting engine measurements is dis- 


cussed. 


iie problem, of computing Opimal reference trajectories is dis- 
cussed in Chapters V and VI. By first deriving the Euler equations with 
the calculus of variations and then demonstrating that the set contains 
a superfluous Euler equation, the linear guidance theory is shown to be 
useiul for computing optimal reference trajectories. The superiluous 


Euler variable is a result of optimizing a state variable. 


a 


To provide a comparison of linear guidance with better known com- 


putation techniques, the steepest ascent formulation of the trajectory 


prebléeimis presentedun Appendixes 


Chapter VIT describes the sample mission that was simulated ona 


digital computer to test the guidance equation as a computational device 


for trajectories. Phe results of the tecmare dicen ca 


Oe Conclusions 


On the basis of the analysis and the subsequent simulation. the 


following conclusions were reached. 


i 


The linear guidance method is applicable to guidance of low- 


thrust vehicles im interplanetary Mieh, 


The rapid convergence of the test procedure proves the mecmas 


to be applicable for computing propellant-optimal trajecionies 
For Col transfers, reserve power for guidance is necessary. 


Mathematical descriptions of many optimization problems 
which contain a state variable in the cost, will produce a 
superfluous Euler variable when the classical calculus of 
variations is used. Deletion of this Euler variable removes 
several difficulties associated with solving optimal control 


Pie@breite: 


Sightings on celestial bodies at discrete intervals may be 
combined with continuous measurements of engine perform- 


ance to estimate the state of the vehicle. 


The uncertainty in state is increased if engine measurements 


are omitted. 


Field-free space analysis of CSI transfers provide a reason- 
able approximation for the results 1n a gravitational field. 
however additional study is needed to ascertain the general 


applicability of the method. 


Ge Contributions of the Investigation 


The items in the report which are believed to be novel are dis- 


cussed in this section. 
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Derivation of the propellant-optimal control law by the method of 
Chapter Ill is thought to be original. The crucial step in the derivation 
is the formation of a reduced adjoint set. This permits inversion of a 
matrix which otherwise would be singular. Other approaches to this 
problem usually require the addition of artificial constraints or weight- 
Pie imatrices which remove the singularity but also change the char- 
ewer oO. the cost function. The method of this thesis preserves the 


cost function and produces the optimal control directly. 


The explanation of the singularity on the basis of a superfluous 
euler variable is thought to be novel. The simplification resulting 
imeomceletion of the extra variable, by using the reduced adjoint set, 
may resolve many difficulties associated with optimal control problems 


Sider than the one in this report. 


The equations which combine the celestial sights and the continuous 
Measurements of engine performance are not believed to have been pre- 


viously derived. 


Pinieally, the method of generating the optimal trajectory 1s believed 
to be simpler than methods previously used. The method results from 
linearizing the state equations but using the complete nonlinear cost 


function. 
8.4 Recommendations for Further Study 


The current research has revealed some areas where additional 
work might produce fruitful results and other areas where information 


is lacking with respect to low-thrust transfers. 


The computational method needs to be expanded to the many body 
problem. Preliminary work does not show any monumental difficulties 
associated with such an effort. The primary problem 1s to include an 


ephemeris in the routine. 


Iti appears possible to use the method as a search routine for finding 
families of trajectories, with the goal of defining optimal launch times 


for low-thrust missions. 


a3 


A study of methods to join different segments of an optimal tra- 
jectory should be interesting. In particular the idea of making only 


-one coordinate change between launch and capture is appealing. 


Research to define rigorously the concept of superfluous Euler 
variables is needed. It is doubtful that the singularities in all optimi- 
zation problems can be removed by the method of this thesis, but per= 


haps some can. Knowledge of the general applicability is needed. 


When data on the reliability of low-thrust power plants and thrusters 
is available, a study of the effects of engine anomalies on the probability 


of mission success should be made. 


A study of VTA guidance using the "critical" plane associated with 
ballistic guidance is needed®?. This concept is difficult to visualize 
in the six dimensional phase space associated with the equations of 
motion. However, the existence of a concept analogous to the "critical 


plane’ idea might provide additional insight to the guidance problem. 
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eee ee 


Trew oie Nee oD =P SPACE 


Peele introduction 


In this appendix the equations for point-to-point transfer of low- 
meee rockets in field-free space (FFs) are presented. The optimal 
transfer problem for a constant power VSI rocket is a straight forward 
me@cation of the calculus of variations. Examples of this transfer are 
mregairequently in the literature. It is reproduced here for compleie- 


mess. 


Foimt-to-point transfer of a Col rocket is more difficult to compute. 
The equations are derived in this appendix. Results of the derivation 
were applied to the sample mission in the thesis to ascertain the validity 
of using FFS predictions in the gravitation field. FFS analysis is found 
to furnish a reasonable approximation for CSI transfers, however the 


analysis is quite tedious. 


The author is indebted to Mr. Neal Carlson for checking the deriva- 


tion and for suggesting different methods of approach. 
fee Constant Power Transfer 


Assume that in FFS we desire to traverse the disiance L in the 
Pome I such that the vehicle begins and ends at rest and such that the 


acceleration integral is a minimum. 


That is 
an 2 
EanGnnetic ies ee ect (A-1) 
0 2 
a 
subject to f vdt=L (A-2) 
0 
and the boundary conditions 
v (O} = 0 
v (T)= 0 tae 
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Forming a functional IF, one obtains 


ie T 
P= | Gira (oes (A-4) 
ae 0 


Direct application of the variational technique produces the optimal 
velocity schedule from which the acceleration is easily derived, 
Ge 5 Zn 
a ae (A-5) 


T T 


Evaluating the acceleration integral, obtain 


2 
ae oe (A-6) 
iy 


From equation (A-5) the optimal initial acceleration is 
a= (A= 


A plot of equation (A-5) is presented in the discussion of Chapter II as 


eine 2 —C. 
A. 3 Constant-Specific-Impulse Transfer 


The well known equations for a conventional! rocket must be inte- 


grated in the analysis of a CSI transfer. That is 
v=cln MR (A-8) 


Since the variational approach is degenerate for this type of optimiza- 
tion problem, ‘we shall simply apply the known result from the Hamil-~ 
tonian approach of Chapter III, that thrust is either full on or full off for 
Csl vehicles in a linear field. Integration of equation (A-8) is Over =i 


initial and final thrusting periods with a coast in between. That is 
ba ae 


L= f v(tyat+ f v(t) dt + v (t,) (t, - t,) (A-9) 
0 ie 
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subject to the boundary conditions (A-3). This rather lengthy compu- 
Paion reduces to 
1 2 
L=c(t, - t,) In |———] tat (A-10) 
Z 1 or 
L=- mt 
1 

where ao in the initial acceleration and m is the normalized flow rate 
merc . a = 1). Some additional algebraic manipulation allows (A-10) 
to be cast in the form 


2 


Palle 
1) 


ic 





(Ary) 


=~ (ar +m -1)1n m, + (1 - m 


where 


m, = 1 - mt, (A-12) 


From the work of Chapter I] observe that for ma positive number 





a =me (A-13 
O 
p= me? (A-14) 
2 
Thus 9 
3 ae 
m = (A-15) 
2p 
c= ae (A-16) 
a 
O 


Equation (A-11) may now be written in terms of initial acceleration, 


power, mass, and the dimensions of the problem, L and T. 





a ie aS qT? \ 3 
5 a ee gir Inm, + (1-m,) 

4p 2p } (A-17) 

in the VSI case ao was a . ler thesCS!) initial acceleration be an 
a 
unknown multiple of that value 
ao = xX a (A-18) 
ap 


on 


Further, define a parameter R: 


2 
R= ea (A-19) 
pr 
Substituting (A-18) and (A-19) into (A-17) one obtains 
a Pe x? = -(R x? + mi - }inm, + (1 -m,)° 
(A-20) 


From the condition that the velocity change must be the same for both 
thrusting intervals if boundary conditions are satisfied, the final mass 


Lamoslrs 


) (A-21) 


Therefore if m, is a Maximum, mass ratio is a Minimum and likewise 


ioe Une vaeeeleroplom Mie no 








om 
Using equation (A-20), solve for x such that Be 0. The result is 
Ox 
om, 
= 0 (A-22) 
Ox 
when 
Rx(5Rx+2Inm,) = 0 (A-23) 


x= O4S 2 iriviel solution and the desired solimicon 1- 


In m (A-24) 


2 in MR (A-25) 
R 


* 
li 


lf equation (A-24) is substituted into (A-20) the result 1s 


a ae 


3K 


et Sm = 0 


)? + (1 ~ m4) In m 1) 


| ] 


(A-26) 


The transcendental equation (A-26) is the relation between m, and R 


for optimal CSI operations in field-free space. The optimum initial 
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acceleration is given by (A-24) or (A-25). Equation (A-26) is plotted 
in Figure 2-d. It is observed that CSI control always uses more pro- 


Sellant than VSI control. Also plotted in Figure 2-d is the CSI curve 


for = fre oo 8) = ins valyevot ao is considerably more expen- 
ae 
sive than a transfer witha _=a . Figure A-a contains plots of 
O O (opt) 


m@iewrequired coast time for CSI control. These are obtained by solving 
equation (A-10) for (t, - t;)/Tin terms of m, and R. Again the work 


involved is substantial. 
Observe that for the case where 


m, = [ere (A-27) 


where € is a small quantity, equation (A~27) may be substituted into 
/ 20) and the resulting expression solved for m, by neglecting higher 


order terms. 


m= 2 (A-28) 
16 
By approximating In m, as 
In ee ne (ely pee) (A-29) 
one obtains 
a 3/4 (A-30) 


Mnereiore, for transfers such that the propellant consumption is 
mall with respect to the total initial mass, the optimum initial acceler- 
ation for CSI vehicles is approximately three-fourths of the correspond- 


ime Optimum initial acceleration for VSI vehicles. 
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APPENDIX B 


ADJOINT SYSTEM OF DIFFERENTIAL EQUATIONS 


ise) Introduction 


mdecepurpose Of this appendix is to review briefly the concept of the 
adjoint system of differential equations and to illustrate its use ina 
Satiple puidance problem. Although the concept is not new, it has not 
been used as a Standard engineering technique in guidance and control 
problems until very recently. For solution of the two-point boundary 
value problem which occurs in the navigation and guidance of space- 


craft however, the method of adjoints is a particularly useful tool. 


The development in this section is not intended to be rigorous and 
exhaustive but merely illustrative of the method used throughout the 
mesic. The material is taken primarily from the lecture notes of 
meoressor Frank D. Faulkner of the United Siates Naval Postgraduate 
School. 


See Wiethod of Adjoints 


Consider the following ordinary scalar differential equation. 


SO ee es a fe) (B-1) 
mare may be used to define the operator 
2 
cc ee eae (B-2) 
dt dt 


where f(t) is an arbitrary but known function. In a manner closely re- 
lated to the method of variation of parameters, form the integral 
a a 
J “a A.L(x) dt = J A f(t) dt (B-3) 


where A is an unspecified function which will be chosen later to satisfy 
certain boundary conditions in addition to a functional relationship. In- 
tegrating equation (B-3) by parts to eliminate the derivative of x from 


the integrand, one obtains 
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7 OT 
J=(Ax+3Ax- Ax l|)+f x(A-3A +2A ) at 
0 


(b=33 
If A. is chosen such that it satisfies 
a 
Ee) = (ea (B-5) 
dt dt 


then the definite integral J is a function only of the first term on the 
right side of equation (B-4). The operation of integrating by parts to 
eliminate the dependent variable from the integrand defines the adjoint 
operator L*. If L* = L the system is said to be self adjoini and has 
some very useful properties. However, these are of no concern at the 


moment. 


A general solution to the adjoint equation (B~5) is 
Ke Coe Seer 


1 9 (B-6) 


Suppose that the desired quantities are x(T) and x(T) and that x(0) and 
x(0) (i.e. two constants of integration associated with the original dif- 


ferential equation) are known. Equation (B-4) may be rewritten as 


ae 
[A+ (3A -A)x}=[AZ+(3A-A)x] + f Af(t) dt 
t=T t=0 0 
(B-7) 
If one specifies that 

IACI) (Bae) 
3 A(T) SAGO (B-9) 

then equation (B-7) gives for x(T) 

| T 
x(T) =[Ax+(3A-A) x] + f Af(t) at (B-10) 
t=0 0 


The constants in equation (B-6) may be evaluated from the boundary 


conditions of the adjoint variable, equations (B-8) and (B-9). Designate 


this solution A ,, By performing the algebra one obtains 


Pe 2 oes (et (B15 


] 
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Equation (B-10) may now be solved explicitly for x(T). 


a 
x(T) = (2e°? -e7') x, +2(e*) -e') xy + f A, f(t) at 
| (B-12) 
To obtain x(T) it is necessary to specify different boundary condi- 


tions for A.. Designate the solution satisfying these conditions as A 9° 


Nepal) 0-220) (B-13) 


3 A.,(T) - IN (ear (B-14) 


2 


By solving equation (B-6) subject to the boundary conditions of (B-13) 
and (B-14) one obtains 
iN2e arial - 2 (t-T) 


2 (B-15) 


Inserting this result into equation (B-7) yields for x(T) 


a 
x(T) = (e7" ~ 7") x, + (2) -e 7") x, + fA, f(t) at 
L (B-16) 
The preceding development illustrates the general technique to be 
used. However, in order to more closely approach the formulations 
used in the thesis, it is convenient to write the basic second order dif- 


ferential equation, (B-1), as two first order equations. 
oe ay, Gas lap 
Va ee ea Lt) (B-18) 


Associate with the two equations two sets of adjoint functions A ‘1 and 
aN where the subscripts 1 and 2 refer to the associated equation and 
the subscript i will refer to the boundary conditions which will eventually 
be assigned to the A's. The integral J may now be formed. 

e at 


i ne N(x -y) + A,, (y + By + 2x) dt “2 Sepeic yact 


(B-19) 
Integrate by parts to obtain 
ie Hla 
T=(AiyxtAiyy) | - x (Ai, - 2s) t y( Nin - 3Aja t+ AG )at 
(By) 
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In order to eliminate the integral term it is clear that the adjoint vari- 


ables must satisfy the relationship 


Aa 2h 8 (B-21) 
furs -3 Ain tAj, = 9 (Boe 


By reducing the two equations to one second order equation and solving 


subject to the boundary conditions 


ea ae (B-23) 

een) = (B-24) 
and 

A, (T) = 0 (B-25) 

Nee (ae (B-26) 


it is apparent that the adjoint system of equations (B-21) and (B-22) is 


identical with that obtained by the first formulation of the problem. 


In general, if any set of first order differential equations mayvaae 


written in the form: 


xy xy f(t) 
2 ee = > (B-27) 
er & sal) 
or equivalently 
x - Ax = f(t) (B-28) 


Then there exists a set of adjoint functions satisfying the relations 


A +AA=O0. (B-29) 
A(T) = 1, (B-30) 
such that 
il 
x(T) = A(t,) x(t,) + J A(t) £(t) at (B-31) 
t 
1 
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where A(t) is the n by n matrix of adjoint functions, Ais then by n 
matrix of coefficients in the original set of differential equations and 


I is the n by nidentity matrix. 


Unfortunately, in most problems of interest the adjoint functions 
Cannot be obtained in closed form and it is necessary to integrate the 
n? equations of (B-29) backwards in time from T to ty numerically, 


using (B-30) as the initial condition. 
fmoeUse of the Adjoint Method in a Simple Guidance Problem 


In this section the set of adjoint functions for a power-limited 
Space vehicle in two dimensional field-free space will be developed 
meeetically. The purpose is twofold: 1) to further illustrate the tech- 
nique of section B. 2, and 2) to provide a simple analytic example for 
use in illustrating the guidance methods expounded in the body of the 


thesis. 


The fundamental guidance equation for thrust-limited vehicles has 
been defined in earlier sections as 
ne 
6s,= A, 6s,+ [ AB 6f dt (B-32) 
=k a ae = 


For an actual interplanetary transfer the A. matrix will be obtained 
pmaiuitanecously with the desired trajectory. In any event, if the tra- 
jectory is known the adjoint functions are also known or may be easily 
determined. To illustrate this assume a transfer in field-free space 
such that motion is along the positive x axis commencing from rest at 


the origin att=0Q0. Att =t, the thrust is reversed to effect rendezvous 


R 
with point Xp at time t = Ups Seer iouire a=. 
7 
t = 0 t = tp t= ty 
—— — 
0 
Pa Tete ee leyiecrel 
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The engine produces constant thrust along the nominal trajectory (which 
is not optimum). Presumably the engine must be controlled in some 
way to provide steering but it is not necessary to specify the method 


of control in order to determine the A. matrix. 


The equations of motion and constraint are: 





Kev, (B-33) 

=v B-34 

. y ( ) 

as 2p) em (B-35) 

eo B-36 

y ( ) 

m= -2p/c? (B-37) 

where ; 
p = —*t~ = (exhaust power) (B-38) 
a 


c = exhaust velocity 


Since the velocity and position along the nominal trajectory are easily 
computed as functions of time by direct integration, they may be con- 
sidered as known functions and we may proceed directly to the varia- 
tional equations which are of primary interest. Taking the variations 


of Equations (B-33) through (B-33) yields 





oe = Oy (B-39) 
x 

by = bv, (B-40) 

by_=- “P 6m + = (2) (B-41) 
cm Im S 

bv = 0 (B-42) 

y 
6m = 0 : 5(2P ) (B-43) 
& 


These may be formed into the matrix equation (B-28) 


bs = A bs + Of (t) (B-44) 
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Or 


ioeet co 0 
oO 0 

és = 0 0 0 0 -2p/em*{ fs+ <¢ = 6(2p/c) 4 (B-45) 
000 0 0 0 
000 0 0 - 8(2p/c7) 


The column vector containing the power and exhaust velocity varia- 
tions corresponds to the arbitrary forcing functions, f(t), of section B. 2 
and to the product Bof in the general guidance problem. They may be 
written in any manner convenient for the investigation at hand. They 
are of no further concern at the moment and will be carried along as 
B of. 

The adjoint functions may now be determined by direct integration 
of equation (B-46). 


Ae = O- (B-46) 


A(T) = 1, (B-47) 


Expanding the above and integrating the twenty-five equations one ob- 


tains 





~ m al 
it if 
[GSE en Sh) o> ag 
m LG m m 
\c me = ah sma bt a 
Mp f 
0 1 O (tpt) O 
1 1 
A= [( - ae ee 
0 O 1 0 ~C ! 
2 1 1 
ee a aie ee 
(a mM, m R 
0 0 0 i 0 
0 O 0 Q 1 


(B-48) 
Thus the adjoint set has been determined for a simple case and may 


be used in the study of guidance techniques. 


nO, 


APPENDIX C 


PROPERTIES OF THE FUNDAMENTAL 
GUIDANCE EQUATION 


Cty ouimiiary, 


The contents of this appendix include discussions relative to inter- 
changing the control vector and relative to properties of the adjoint set. 
In the absence of control perturbations the adjoint set alone describes 
the effect of state perturbations and thus may be called the "state tran- 


sistion matrix’. Several of its interesting properties are discussed. 


Proof of singularities in the guidance equation is presented in 


section C. 5. 
C.2 General Remarks 


To provide a better understanding of the methods used in t(hese@nm 
tion of the equations of motion, and to show the relationship between 
trajectory determination and guidance it is beneficial to examine vari- 
ous formulations of the variational equations of motion and the con- 


straining equations. 


A by-product of using the adjoint method with the calculus of varia- 
tions to find an optimum trajectory is the solution of the adjoint set of 
equations. Physically the adjoint set represents ''sensitivity' coeffi- 
cients or ‘influence'' functions which permit the investigator searching 
for a trajectory to determine the changes in trajectory parameters 
which will move the solution in the direction of the desired optimum, 

It is interesting to note that along the optimum trajectory the Lagrange 
multipliers are linear combinations of the adjoint variables. From 
the viewpoint of one searching for a trajectory, the adjoint set has 
served its purpose once the optimum trajectory is determined. How. 
ever, for the guidance analyst the adjoint functions serve a most useful 
purpose by showing the effect of spacecraft perturbations on the final 


state of the vehicle. 
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feo interchanging the Conirol Vector 


The equations of motion and the mass rate equations, expressed 


in a nonrotating frame with origin ai the central body, are 


ia (C-1) 
v=-Erta (C-2) 
3 
lie 
m= g_ (m, a) (C-3) 
or m=g_ (£) (C-4) 


The variables m and a have a functional dependence which must be ex- 
plicitly taken into account. In the subsequent discussion only the VSI 


case will be derived since CSI control follows an analogous argument, 


The relationship between m and a may be written in any manner 
which is convenient to the argument at hand. If one is interested in the 


effect of acceleration changes irrespective of their cause then write 





(C-9) 
2p 

The analysis of Chapter I] shows that equation (C-5) 1s the correct 

form for mass rate. By taking the first variations of equation (C-1) 


through (C-3) the state variational equation is obtained 


Be See aes a (C-6) 
where 
on nO 
aC ROO 
a 3 9 
Ce On c7 





0 
Os 
ae I, (C-8) 
-m* sl 
a 
_p - 


where the pre-subscript denotes that ais the control. The submatrices 
of A and B are the three by three null matrix O23; the three Dy iia. 
identity I, the three by three symmetric matrix of gravitational grad- 
ients G; and the three by one null vector O. If the adjoint set Torjeqmes 


on (C=6) 1s formed such tiat 
A=-/NA (Gag) 
ene (t.) = | (C-10) 


thenrthe Solition. = must pave We 1)6em 


ae ee oe (Cree 
ae ale 
oT of Ag, 


where the matrix is partitioned into three by three submatrices, three 


by one vectors and a scalar,/\ 49. 


The form becomes apparent when equation (C-9) is expanded and an 
initial condition is applied. Notice also that for the form (C-I11) te 


exist it is not critical that (t,) = I, only that at some time t,/U(Gjmm 


The expansion of (6-9) yields 


Nae AG (C-12) 
Ay? Ay Ce 
re (C1 
Kgy hag (Cay 
Nop "Soy (Co 
eee (Camm 
Aah c-9 
Ag2> Agi (Ca 
Ne : a2 Lae (C2 





——e 


The form of equation (C-11) is assured if at some time /., = I. 


Now rewrite the functional dependence of a and m as 


f 
a, ie 21) 
(2 
ne (C€=22) 
2p 


where again f, m and p are per unit initial mass. The A and B matrix 


must be modified to 


_— 


Our <1e 
a 
A= [-G O, -= (G=25) 
eons 
Os 
2 es ; 
> | als ce 
ae 
“Pp. 


where the subscript denotes that f 1s the control. The state variational 
differential equation is 


8s = EOS 1: 


; Bot (C-25) 


f 


Clearly the optimal trajectory must be independent of the particular 
control variation used in the computation scheme. That is, if the varia- 
tional problem for computing trajectories is solved using f as the con- 
trol, it must result in the same trajectory that would be obtained if a 

is used as the control. Thus, two different formulations are available 
mime may be used at the discretion of the investigator. Since the cost 
has been specified as aye it is more convenient to use a as the control 
for computing trajectories. Nevertheless, the use of f as a control has 


important significance. 


The adjoint set which corresponds to using f as the control will now 
we ceierimined, Mxpandine (C—O sine (C-235) One obtains difierences 
from the set (C-12) through (C-20) only in the last column of the solu- 
men /\, lnat is 





a 
® — 
ne ele (C-28 
‘ = 
A. 93 “29m (Co 
NN. 33 = 0) (C-28) 
For boundary conditions /\. (t.) = I, the solution has the form 
hoe ee 
t= | Noi A223 (C=29) 
OF Fou 


The solutions a and p/\ differ only in the last column. The renaauaime 


numbers of the array are identical for each point in time. 


The physical significance of this difference is apparent from the 
following physical reasoning. Consider the adjoint set as an array of 


partial derivatives. IfJA. (te) = I then the last column is 








Ox (te) 

@ m (t) 
a dv (ty) oa 
D 23 3m (t) | 
A\.33 om (t, 

dm (t) 


If ais the control then one should expect N13 and /\. 9 to be null veges 
tors, Since for a given acceleration program variations in mass can- 
not affect the end point in position and velocity. However the end point 
for mass must change because it requires a different amount of pro- 


pellant tee dilierent size velnicics, 


If f is the control and the thrust program is given, then changes in 
mass will affect the end points in position and velocity due to the change 
in acceleration. For a given thrust program the final mass, however, 


can change only if the initial mass changes, therefore /L 33 7 ile. 


1A 


It is concluded that the control vectors 6f and 6a can be interchanged 
inthe guidance equation simply by changing the last column of the ad- 
joint set and changing the matrix B. Along the optimal trajectory all 
of these functions are known thus the interchange can be easily made at 


e@esoiscretion of the investigator. 


The purpose of interchanging control vectors is simply stated. Com- 
putation of trajectories is most conveniently carried out with a as the 
control, due to the formulation of the cost. Assume however, that it 
fewcesired to investigate the effects of changes in engine performance 
or to investigate variables other than a as quantities to be measured. 


men it is convenient io write, from Chapter I] 








pie rina: (Cas) 
then 6f = méc+c 6m (C-32) 
6m | 
or of = [c m1I,] (C-33) 
f=[ce - 
Further, let f£=-=P ¢ (C-34) 
Cc 
ae 20 2ect 
then Beso yan Te, Pea oe dl) 1G (C-35) 
7 ce cf 2 of S 
or 
9 2cc 
Ji) geen Boomer bc oe) 
< eS oa 


If in equation (C-25) or in the guidance equation, of is replaced by 
ener {(C-33) or (C-36) it is possible to study the effects of the par- 


ticular variations represented without computing a new adjoint set. 


The ability to interchange control vectors at will, without requir- 
ing a complete new set of computations makes the fundamental guidance 


equation a very powerful tool. 
C.4 Properties of the State Transition Matrix 


If the fundamental guidance equation is written for two different 


times t = t and t = ue and subtracted, then 


Uj 


A(t) 68, = A(t.) 68+ ! AB 6a dt (C=3m 

al t. 
bs, = Tj, 6s, + A’ (t,) [a2 ba at (C-38) 

where | 
T= AS (ty) A(t) (C-39) 


In the absence of control perturbations, Tr completely describes tie 
transformation of the state vector at t. into the state vector at t.. Thus 
T.. is called the state transition matrix’. It is necessary to show 
that a is ee It is desirable to show that it canbe inverted by 
inspection and that Te = att This last property is readily observed by 
taking the inverse of equation (C-39). 

-1 


-~1] _ -] 
to eth Ae (C-40) 


i) = ee ee (C41) 


1j ia "i 1 


To examine the other properties, it is necessary to look first at/.. 
Prom section ©. 3 wriley\ in its seneralaonm 
Nyy Aye Ais 

A= 1No1 Ago Ags (Cra 
OF OF N33 

Equation (C-42) is to be understood as representing both the adjoint set 


with fas control and the set for a as control. The last column takes on 


the value appropriate to the control. 


Now, partition /\ and consider only that portion defined by 


Sg eo 
Hi = (C-43) 
No, Sa 





Mimere the Subscript refers to time t.. The differential equations for 


mis Subset of /\ are reproduced from section C. 3. 


2 shee Sale 2) 
Nit Ady ee 
on “/\. 92% (ote) 
N97 "N94 | (Sea) 


Observe that these equations are independent of the control and are 


completely specified by the physical path and the boundary conditions. 


oie O-13s and C-16, et may be rewritten as 


~/\ vy 
H+ “12 12 nae 
Ngo Noe 
Now define a new matrix en 
Nig Nag 
Eee ows . (C-45) 
J Ni2 N22 
a and el abe mellated by 
a HOP (C-46) 
j j 
where 
O -] 
Pp = 3 3 (e227) 
I Og 


oy (C-48) 
eee (C-49) 


pince in equation (@—4) t is an arbitrary time along the trajecton,, 
the subscript may be dropped and ae considered as a time varying 


Matrix. 


3 


He > Serer (C-50) 


DiMerentiating: 
H*+ = HP (C-51) 


From equations (C-12) through (C-16) and (C-44) it is possibiemg 


write 
H + HA* = O (C-52) 
where 
oO, I 
ae=| 2 8 (C-53) 
2 O. 


substituting (C-52) into (C-351) and using (C-50) one obtains 


af ae 


H*+ = +H*" PA* P (C-54) 


Multiplying out the known product yields 
parle heen 
PAs Pak (C-Siaq 
since Gis a symmetric matrix. Thus the transpose of equation (C-94) 
1s 
H* _ A* H* (C-56) 
Comparing (C-52) and (C-56) it is observed that H and H* satisfy the 


adjoint-fundamental relationship mentioned in Chapter III. That is 


d 


an (H H*) = O¢ (C-a%@ 
or integrating 
io eek C-58 
f ee ( ) 
where 
lal (te) = He S I, (C-59) 
Substituting equation (C-50) into (C-58) and using (C-49) obtain 
-PH, = - H, one (C-60) 





misime (C-959) 
H, PH. = p (C-61) 
iy inecrix Hi exhibiting the property of equation (C-61) is said to 
be symplectic and its inverse is given by 
lap os isles (C-62) 
The inverse of the submatrix H of A. has been obtained by rearrange- 
ment of the components. To obtain A two cases must be considered: 


iieeadjoint sets for the control a and for the control tf. 


If A. = aed Ns then 
H Of, 
ae (C-63) 
ee 
where QO, ise omeOmpOnens Null vectoreu fF urtiner 
1 [HQ 
a) oes (C-64) 
Ore t= 
BU Us 
fi /\. = p/\; then 
= A ¢ | 
Pe Ga (C-65) 
a 
CO; 
A143 
where na denotes a= 
a A293 


Equation (C-65) may be written as the product of two matrices. 


i - o a6 - ~6 (C-66) 
oF 1 OF 1 
Then ea Og -] N¢ 
i . (C67) 
Ony 1 naan 





Ns (C-68) 
O% 


since the second term in (C-67) is its own inverse. Therefore _/ has 


now been inverted for all cases. Now consider T;.. 


J 
ys ae (C-39) 
1) - J 
For A = au , use (C-63) and (C-64) 
eae 
i te 
tT. 8 [2S (C-69) 
a ae | /\ 33; 
8 | N33 i 
For/. = pits use (C-65) and (C-68) 
Lee en 
Hy Hy 1-H (Ai A 6; 
Te tee = Sain eee (C-70) 
Q5 1 7 


Finally it is necessary to show that ar is nonsingular. To do this 
it is sufficient to show that the determinant of i: et is never zero and 


that /. 39 is never zero. 


From equation (C-20) it is apparent that if NN. 33 (0) = 1 it increases 


monotonically, thus is never zero. 


Now examine the determinant of 
Eien 
Prom (G—0 1) 
jepie Siden Pe (C-71) 
Since the determinant of a product is the product of the determinants 


(Geno (C-72) 
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from equation (C-62) 
1 ae 


det  =det{-PH P) (C-73) 
aise = 6h (C-74) 
since det (-P) = (-1)" det P = 1 (C-75) 


where k = 6 is the order of P. 


emilee the Subscripts 1 and ] are completely arbitrary and since H 
is continuous, 


1 1 


det H, Hi, =detH H=1 (C-76) 


Thus ar is nonsingular. 


©. 5 Singularities of the Guidance Equation 


Solution of the guidance equation for a control program which will 
satisfy certain terminal conditions invariably results in the inversion 
Sieeimiatrix integral expression. In section 3.6 the solution for FTA 


guidance requires such an inversion, namely 
tf a, 
mt- |f AxBB i A*?at (Cana 
oi 
In this section the proof will be given that M is not singular for 
(tp - t,) > 0 but that the corresponding matrix denoted by M_, with/A. 
imeplace of /\*, is singular in the case of unconstrained control pro- 
grams, for all values of ty along the optimal trajectory. We shall pro- 
Seeq tO prove that M, is Singular and by induction show that M is not 


Singular. 


ime tOllowines proot of the Singularity in M, was suggested by Dr. 


wames Potter. 
If Mo is singular then there is a nonzero vector Py such that 


Nor — 0 (C-78) 
Oo =o 


ie 


M, is symmetric thus it is also true that 


Po M, = 0 (C=3ay 
and 
oi M p = 0 (C-80) 
=O Oo 
This implies that 
Ue 
0 = p. ABB A’ pdt (C-81) 
i 
1 


However equation (C-81) is the integral of the square of the length of a 


vector 


Hi | Be Ni p | Gi (C=875 


Thus for Mo to be singular it is necessary and sufficient that the vec- 


LOG INte prance OL (C762) beaZere: 


B’ (t) A’ (t)p, = 0 (C-83) 


Tor all oi i te. 


ia veCcror Po exists for the optimum trajectory such thatveqtas 
tion (C-83) is satisfied may be shown by application of the Pontryagin 
maximum principle for unconstrained control. The Pontryagin maxi 
mum principle asserts that if the final value, 5, of some combiftanias 
of the state variables, x, are to be a maximum or minimum with re- 
spect to the control variables, y; that is, 
minS=d" x, (C-84) 


i 
then the desired path will be one such that the gradient of the state 


vel@cily.. x4 inl the veclor space er -allacmissi ole. y, 


i g (x, y, t) (C-85) 
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will be orthogonal to a vector p which satisfies the differential equa- 


tions adjoint to x = g (x, y, t) and such that the final value of p is 


pes av, (C-86) 
where d. = 0 for components of the state vector which are fixed at Ue 
and a O for components of the state vector to be minimized or max- 

ft 
mnized., 


Thus to satisfy equation (C-84), it is necessary and sufficient that 


T 0(g (x, y, t)) 


p ee = (C-87) 

oy 

where 
ae gs (x, os t) (C-85) 
8 (g(x, y. t))] 

oo p (C-88) 

oes 
Prpz-dty, (C-86) 


A proof of the maximum principle is given in Appendix D. 


Application of equations (C-86), (C-87), (C-88) to the guidance 


problem in this thesis reveals that 


o y V 

pee - ee (C-89) 
1 Ong ey 

p=-A’ p (E290) 

p' B= 0 (C-91) 


where ee and B are the matrices of partial derivatives in (C-3938) and 


We-o/) respectively. 


Since p satisfies equation (C-90) and also At satisfies (C-90), a 


solution of pis 


ee ea, (C-92) 


ane 


Then substituting into equation (C-91) 
ie 2 


Pours lt t,< fot 


ty is arbitrary 


f Thus Pe is a null vector of M |. Furthermore, since 


M.(t,) Pp = 0 (C205) (C-94) 


To show that M is not singular it is sufficient to show that every 
MULE Ve Clo Ol M, satisfies the boundary condition on p, equation (C-86). 
This sufficiency condition is easily verified by noting that the final 
eulerior signe Oi De is invariant and that M may be obtained from M by dee 
leting the last row and last column of M In general, if /\ is inter- 
preted as the nonsingular transformation that transforms ) into p 
and if p! is a null vector of B because of the component Pip? then deleuws 
ing the ith row of A. to produce /\* destroys the transformation which 


carries Pig AAC 


The fact that only the last row of /\ need be deleted is a con- 
sequence of the fact that only the last component of the final state vec-— 


tor, m,, has been maximized in obtaining the optimum path. ThereiGimes 


fe 
we assert that if every null vector of M, satisfies equation (C-89), 


the M obtained in equation (C-77) is nonsingular. 


Let us prove then that all null vectors of M, satisfy the boundary 
condition. For the problem in this thesis B is a seven by three miavma 
of rank three and contains zeros in the first three rows. By virtue of 
the rank, the three columns of B are linearly independent. Furtiew 
since Bis of rank three there are four linearly independent vectam. 
which satisfy (C-91) and thus are null vectors of B. Three of the four 
vectors may be chosen to satisfy (C-91) by virtue of the three rows of 
zeros in B. These vectors do not satisfy (C-92) and (C-93) and there- 
fore are not null vectors of M |. The fourth vector must be a null vec- 
LOreot M, Since the three columns of B and the four null vectors of B 
form a basis in seven dimensional vector space and six of the vectors 
are not null vectors Of Mo: This vector is unique and must satisfy 


(G-o7)} throug ( @-93)), 
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The proof is thus complete since all vectors may be written as 


linear combinations of the vectors in the basis. 
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foe PEN pee 
PONTRYAGIN'S MAXIMUM PRINCIPLE 


D.1 Introduction 


In this appendix a derivation of Pontryagin's maximum principle 


is given which is applicable to the problem of the thesis. 


Consider that the cost function S is some scalar function of the 


final state variable, x It is always possible to cast a problem into 


Chaise lori Oy. a. variables. The control vector will be denoted 
by a generalized vector u, which is bounded by certain constraints, mm 
addition, the solution must satisfy boundary conditions on some of the 
state variables. The desired solution is the optimal control, ue which 


will maximize the cost and satisfy all constraining conditions. 
| Da Wi ot) Fe aheoybbecwe emanhetonholk= 
Assume that the state variables may be written in the form 
x= g (x, u) (D-1) 
and that the solution is to satisfy the final boundary conditions 


$ (x-) = 0 (D-2) 


The cost function may be written as 


ae 


_ ar 
S=d  x,-tyv @ (D-3) 


where dis a known vector and Y is an unknown constant vector. 


Then small changes in the cost due to changes in state are given by 





a9 
tae Oca ce (D-4) 
Se - 3, x 
— f 
where 
GLO” 
two 0 (D-5) 
OX. 
if 
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memation (D-5) is 1 for fixed boundaries; 0 for unspecified boundaries. 
mene variables adjoint to (D-1) are given by 


: le 
oe ee (D-6) 
Je ee \ } We 


Cc 


@ 
[x 


I'S 


(t-) = De (D-7) 


iiiesvector p is frequently called the ‘costate’. 
mosume that an optimal control, u’, is known. Then for small 


changes 6u, the state must change by 6x. Then 


6x = g (x° + 6x, u? + bu) - g (x, u°) (D-8) 


loa 


Consider the term 


Ox (D-9) 


Applying (D-8) one obtains 


p’ 8x =-2- (p" 8x) - p' [g (x° + 8x, u° + du)-g (x°, u)] 
o (D-10) 
Integrating (D-10) yields 
Pe. tee 
Jp oxdt= pox | - f pe tg (x° + x, u? + bu) - g (x®, u®)Jat 
0 0 0 
(D-11) 


If the initial state boundary is fixed, 6x = 0 and the first right hand 


term is 


Taos 5x (D-12) 


provided Pr is chosen such that 





20) 
pe (Dee) 


p 


OX ¢ 


The right side integrand may be expanded in a Taylor series around the 


optimal trajectory if the variations due to control and due to state are 
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Se€pardiec. == har is 


g(x + 6x uF 6U) 2 Gem gee 


: u~ + du) - g e. u) 


4s 2g (x°, u> OU) [Oe er re Maa itiele te 


a (D-14) 


lf the adjoint relationship, (D-6), is used on the left hand integrand 


oi (D-11), the term may be rewritten as 


Uf 


ee a f Be tg (x°, u-) ] 8x dt (D=a) 
O° 0x 


Then substituting (D-12), (D-14) and (D-15) into (D-11), the resulta 


i 
it 
6S= f p> [g (x°, u° + bu) - g (x°, uy] dt 


0 bs = 


rs otf g [g (x°, uv + bu) - g (x°, le ox dt (D-16) 
lax a | 

Consider the first integrand in (D-16). If Sis a maximum, then any 

allowable change 65u must cause 6S to be negative, Therefore for all 

as te 


p- [g (x°, u? + du) - g (x°, u 20 (D-17) 


In considering the second integrand it is argued that admissible 
control changes 6u, can produce only a small variation in the deriva- 
tive term such that its product with 6x is second order or higher. Con- 
sequently the second integral and all higher order remainder terms of 


(ete) aberne slectee: 


Therefore (D-17) represents the sufficient condition for an Opti 


A necessary condition is 


p [g (x°, uw? + bu) - g (x°, u)] < 0 (D- 18) 


The implication of (D-18) is that for all points on the optimal tra- 
jectory and for all admissible functions u, pg must be a maximum. 
Replacing g with x, notice that the scalar product of p and the state 


velocity must be a maximum 
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ee ce 


If the control is unconstrained, then a necessary condition that 


pig iS a Maximum is that 
ae 
0 T T & 
a Og Oe (D-19) 
du du 
From previous work 
ae 
Bz— (D-20) 
au 
Thus a necessary condition for a maximum is 
p'B=O (D-21) 


It is possible to derive the Hamiltonian formulation used in Chapter 
III by introducing a slightly different cost function S'. In place of (D-3) 


let 
tf 
s)= f hn(u at+vl (D-22) 
; u uy 


feere (uu) is some function of the control. Then 





tr ag 
Se ee ea ooo (D-23) 
QO Ou Ms Tb bis i 
_ —f 
where again 
oo, 
= rere (D-5) 
Oe, 
if 


Equations (D-6) through (D-11) hold and (D-23) may be written as 


tf 
6s'= f [h (u° + du) - h (u®)] dt + De 5x, (D-24) 
0 
provided that De is chosen such that 
09 
pray (D-25) 
OX ¢ 


Ii D-24) is substituted into (D-11) alongs with (D-14) and (D-lo) sara 
result is 
ee x \ 
6s'= f ¢ p- [g (x°, u° + du) -g (x°, u®)]+[h (u? + bu) -h cor dt 
O 


7 Che aoe oa ihcraiiac (D-26) 


Equation (D-26) is analagous to (D-16). A sufficient condition that $' 
be a maximum is that for all admissible control changes, 6u 


O 


p” (g (x, u° + 6u) - g (x®, u® 


)] + [h (u° + 6u) -h (u°)] < 0 (Da27} 


The implication of (D-27) is that for all points on the optimal trajectory 
and for all admissible functions ie (pg + h(u) ) must be a maximum. 


Using the VSI vehicle as an example 


2 
na) (D-28) 
Z 
Boemausine (M-1) define 
a” te 
H=—+p x (DeZiey 
2 


For the linearized guidance problem of Chapter III, the vector & 
was interpreted as the velocity of the state error at the final time. 
Since a is a state velocity and Ps is the final value of the costate, using 
(D-25) 


a’ eee 
H=—+yp & (B=a0) 
2 
where for fixed boundary conditions 
Me 
a i (D-31 


0§ 
Equation (D-30) is a linearized version of the general Hamiltonian 


formulation, but with a nonlinear cost. 
D. 3 Optimization Criterion 


It was briefly mentioned in Chapter V that the optimality criterion 


for Pontryagin's maximum principle is more workable than the calculus 
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of variations criterion. This remark deserves further discussion. 
From the previous section it may be noted that the criterion for optim- 
ality may be stated: ''For admissible control variations, the change in 
eeeecost Iunction must be zero or negative (positive) if the cost is max- 
imum (minimum). Admissible control variations are those which 


violate neither the boundary conditions nor any control constraint. "' 


In conirast to this, the calculus of variations criterion for optim- 
Meme is derived from the integral of the functional F used in Chapter V. 
Peterring to equation (5-18), for example, the optimality criterion may 
be stated: ‘For admissible control variations and variations in state, 
the boundary conditions must be satisfied and the cost must not change 


to first order. "' 


since the integral of the functional equals the cost, the form of 


equation (5-18), rewritten here as (D-32), 


te otf 
Mich 10) = + f ( yér+( ) bv + 
or 4g, = / 
( )dm+(_ +) ba | at (D-32) 


gives the impression that the criterion must hold for all arbitrary 
variations in state and control. The implication is not true for the 
problem of this thesis and is misleading at best because (D-32) does 
not represent the first variation of the cost. It will now be shown 
that although (D-32) is satisfied for all variations in state and Conroe, 


the optimality criterion can not be satisfied for all arbitrary state 


wariations, 


Let the guidance equation (3 27) be partitioned such that it con- 


sists of two equations 


tf 
for = nN, bs, + f A’Béa (D283) 
ow f jb 
te 
: oe A 7 
6m , = Ay os,+ J A. 7 Béa dt (D-34) 
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Using the expansion of the matrix A_, (D-33) and (D-34) are 


tr © 
“ Det eee ibe Se a + f A'B 6a at (D-35) 
eM f\a1 A22 |, ens 


“ 
Z 
= i m & 3 
6m | = Ngo om | i Az (=) a’ 6a dt (D-36) 
t p 
Using the solution N22 = m*/m* from equation (5-25), (D-36) miayaige 


reduced further. 


ih 


ee 
5m . = a ese f atbeaa . (D=3) 
p 


iL 


For the criterion to be satisfied it is necessary that the final state varia 
5 
tion 6s,vanish. First Siew -\ to take on a nonzero value but let 6m 


SeeleZze corn 1 ene it) 5.5) and (D-37) become 
tf 
O=-&+ f A'Béadt (D-38) 
; | 
. m*, tr 
0 = 0o-—1t f[- al baat , (DEaen 
oe 
From Chapter III it is known that all propellant optimal control varia- 
tions which satisfy the boundary conditions are 
a=BiA * M" (n-&)-a (D-40) 
substituting this solution into (D-39) and expanding yields 
AZ 
C= ae [oP ty-ay-nT arts (D-41) 
p 
However, along an optimal trajectory 
n' M!n-25 = 0 (D-42) 
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Thus to satisfy (D-39) it is necessary that 

O=-n ME (D-43) 
(D-43) does not hold in general. It holds only when the error & is such 
that a control change 6a orthogonal to a can satisfy boundary conditions. 


Now consider a variation in OM ¢ but no variation in position or 


velocity. Then it is required that 


i 
o= f A*B a dt (D-44) 
i 
2 
Qpagunt Poni ag ) 6A (D-45) 
m, amo. 


A general variation 6a may always be written as 


a= Bl Ax! a-a (D-46) 


Inserting (D-46) into (D-45) and (D-44) and solving, one obtains from 
(D-44) 
-1 


Tai Wey) (D-47) 
and from (D-45) 
6m 
ie (D-48) 
m 
t 
Eliminating 7 
dm 
i eee (D-49) 
a 


Therefore there are no variations 6m, which permit the boundary con- 


i. 
ditions to be satisfied and which do not change the cost. 

It is concluded that the calculus of variation criterion and formula- 
tion can be misleading if not treated with much care. The Pontryagin 
principle avoids the difficulty by making no explicit requirements on 


state variations. 
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APPENDIX E 
STEEPEST ASCENT FORMULATION POR GuiDANce 


i l ineroducrion 


In this appendix the guidance problem will be formulated as it would 
appear in a ‘steepest ascent'' trajectory computation procedure. The 
steepest ascent method is characterized by a rigid control covera@iae 
step size between iterations. This step control is to prevent thespmes 
cedure from violating linearity assumptions as it approaches the opti- 
mum. sophisticated techniques have been devised for automatic selec- 
tion of step size. a However, at some point the investigator must use 
his experience and judgement to select a constant or matrix of coms 
stants to insert in the selection procedure for step size. The author 
is convinced that the additional complexity of rigid step control cannot 


always be justified. 
E.2 Steepest Ascent Solution 


For simplicity in illustrating the steepest ascent method, comsig@es 
only the trajectory problem for the unrestricted VSI mode of control. 
In this problem it is desired that the optimal acceleration integral be 


eae eg Ute 


dt (Haas) 





It is further desired that the change in the acceleration progpram™anaas 


one iteration to the next be such that 
2 f : 
kK = f dal at (E-2) 
0 << 


where ke is a constant to be selected at the discretion of the investigator. 


In order to satisfy, to first order, the terminal constraint &, 6a 


must satisfy the relation 


Loz 





te 
-&= f A Béadt (E-3) 
0 


where all variables are defined exactly as in Chapter III. 


Define the new cost function J' such that 
aE 
tf Bo 0 


7 mt 
gz f a+ (a? - f sa" badt)+ n> (s+ f N'B 6a at) 
0 Z 2 0 0 
ea 


where 7, and 7, are constant Lagrange multipliers, a scalar anda 


i 2 

memeror respectively. 
For arbitrary variations in 6a, the first variation of J' must vanish. 
Thus for 


eek oa (E-5) 


where ais a nonoptimal program, 


O = ta ba)! = 7) bal ts a> As 6 (Oa) (i: -6) 


Or 


i el 
- oe ALE ie :, 
Bie Se eee (E-7) 


(] - 74) 


The multipliers may be determined by satisfying the constraint equa- 


Meme. First solve for 7, through equation (E-3). 


2 


a,=M1{(1-m,) &- 71 (E-8) 
where 

ie : 
M= f A*BB’ A® at (E-9) 

0 

tp 

n= | A° Ba at (E-10) 
iG 
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and 














se L - uf 2 
ba=--B' AN wt fe | (E-11) 
(1-7,) | (1-7) 
Sites 2 tOrcolvc dels one obtains 
(1-7,) 
eo 
(l-7,)° = —,——>,_— (E-12) 
1 Ke 2 etm va 
Oi 
ee (E-13) 
(1-7, ) 





where J is the acceleration integral for the nonoptimum program. 





a ) is the step size control which must be deter-= 
1 
mined by judicious selection of a Observe, however, the result of 


using (H-11) to evaluate the integral of al 6a, This term is the Cyoee 


The factor ( 


PeOGuct il ine Cxpancsien Ol Waele by equation (E-5). 


£-—— (25-9'M * 9) (E-14) 
0 l-7 





i 


As the program a approaches a, the difference, 6a, must vanish. 
Likewise the error must vanish. From (E-14) either the step size 
must approach zero or the term (2J - ni M7 'n) must approach Zeresee 


both. It is easy to show that both must approach zero. Consider 


=e 6 Vel eere (E-15) 


|e 


Then 


Uf 
fatadt=23-y'M'y (E~-16) 
a2 3 0 0 
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elearly 


f ataadt>o (E-17) 
. 2 ee 


with equality holding only when 


1 
1 


a=B? axlu ly (E-18) 


ouation (-18) holds only along the optimal trajectory. Therefore from 


(=13) it is clear that ae ) to remain finite, ie - a M7? £) must 
i 


approach zero at leasi as fast as (2J - n M7" n). since | é | approaches 


zero thisimplies that Pe Mitstaoe medlced TOeZero aS me Optimal patius 





approached. One concludes that "steepest ascent’, as formulated using 


the usual technique ae inherently converges slowly. 


The approach in this thesis is to assume that as Je | paoeaneS small 
the step size will automatically become small without the arbitrary con- 
straint imposed by i Thus 7, is set equal to zero and the procedure 
allowed to converge as rapidly as possible. The results appear to 


justify this procedure. 
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APPENDIX F 
COMPUTATIONAL COORDINATES 


Rl lnteoduchion 


In this appendix the computational coordinate system is derived 
from the ephemeris data of the launch planet and target planet. In the 
chapters describing the guidance equation and its uses, generalized 
vector notation is used and the problem of coordinates does not arise. 


The motion is assumed to be described in a nonrotating frame. 


For computations, however, it is necessary to be more specific. 
The computational frame used in computer tests is a heliocentric frame 
defined by the transfer plane and the initial point. The transfer plane 
is the plane which contains the initial point, the final point and the sun. 
The x axis of the system passes through the initial point, the z axis is 


hortnerly- and the y axic completes ihe triad: 


The objective is to describe the transfer plane and its coordinage 
system in terms of ephemeris data. This data may be given in either 
ecliptic or equatorial coordinates. The description of these systems 
and the transformation between them is available in most basic celes- 


Wale chanics texts and will mot be repr oduced nee 
F, 2 Computational Coordinates 


Designate the unit vectors in any system as i, j, and k and attach 
subscripts to denote the system. For the computational system use 
subscript c. Then, if the launch and target points are denoted by sub- 


Seip = and TP respectively. 


r 
. a Z 

Lee (F-1) 

IL 

Ieee eae 
age ees (F-2) 
~S lr Xx 

ee 
elo see (F-3) 
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Mie transfer angle, 0, in the transfer plane is determined by 


bo 
Ke sin Q = oo ee (i’-4) 
prier cut 
/k sin | = sin @ (F-5) 
The coordinates of the launch (initial) point are 
1s se 
ie (B56) 
Zy O 
c 
The coordinates of the target point are 
Xen rip COS 0 
Yr = Cop sin @ (F'-7) 
<7 G Cc 


Two additional parameters of interest are the inclination of the 
transfer plane to the planes of the launch planet and the target planet. 
If the northerly normals to these planes are denoted by ky and Kin 
respectively, and the angles by ap and Op Pespectively, then 


cosa, = Kr La (F-8) 


k 


eg een ee 


c 


cos a 
where 


Cee ore ° 
k. = ~L (b+ 90°) (F-10) 


le XX, 4 90°) | 


Sap Sa ene 


oie 


(F-11) 
—(T + 90°) | 


To r 
—'T 


ar 


From the preceding equations, numerical values of the coordinate 
transformations (i.e. the orthogonal transformation matrix) between an 
ephemeris tabulation and the computational coordinates can be com- 


puted when the launch date and transfer time are specified. 
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APPENDIX G 
COMPUTER PROGRAM 


G. 1 Introduction 


The FORTRAN program used to test the guidance theory is repro- 
duced in this appendix. The input data format is given as well as the 


units for input data. 


There are a few areas in the program which can be made more 
efficient, since testing is completed, however it operates quite satis- 
factorily as presently written. The author is indebted to Mr. Krupp 


for his magnificent efforts in writing this program. 
eee input Format and Units 


hea@e data cards are required for each problem LO be COMmulLed, 
Mmaestormat for each is 3F20.9. The following sequence and units are 


feecuired: 
ies Initial position (A. u. ) 
2.) Initial velocity (A. u. /day) 
payee Largetl position (A: ul) 
4.) Target velocity (A. u. /day) 


oy pend Gs) yEsiimated y vector. For Vsl control, both cards are 
DUM VeCclome mH Or €5! Control, card muumber 5 1s the null 
VMeclt@reanGdce acd oO should Contain a number of the order o: 


maximum initial acceleration. The units are (A. u. day”). 


7.) a. Number of iterations desired 
b, Maximum initial acceleration in units of 10°*g.. (Bexaimaieiics 
I, 2). ee or Vol problems use any larce number, 
c. Gravitational constant of the central body in units of 
(A. een For the sun this value is 0. 000295912. 


8.) a. Maximum time step desired (days). 


one 


b. Earth gravity (A.u. /day*). This value is 0. 494840, 
Flight time (days) . 





9.) a. Exhaust power ol ecnoine per unitimacc le \* /day?. 
The conversion is: divide power in kw/kg by 3. 4653 X 107. 


Deo W RUG mtd ol Uilad io ete es ue len 0) |e 


The routine will accept as many problems of one type as desired, Ger 
and VSI problems cannot be run together. Submit a complete set of | 


Gata Cards 1orreeaem problem: 
Gas ne POR Mh ane Pooram 


Except for subroutine DERV, the programs for CSI and VSI control 
are identical. There is a subroutine DERV for each mode of ensime 
control. Select the routine appropriate to the problem and omit the 
other. Selection cannot be made automatically with the current pro 


gram. 


The FORTRAN source program is reproduced on the following 


pages. The generalized flow chart is presented in Figure G-a. Figure 


—<.- << 


G-b is the storage map for computer variables. 
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Fig. G-b. Computer storage map. 
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i a a al a 


ASE 

eS 

MAIN 

DIMENSION W( 12s 20)5 Tl7s 8s 60) 

Peer is ([(¥ (ls J)s J=l195s 20)5 [=ls 7s 2)5 W(B» 17) Was 15) 
INPUT DATAs (SEE FORMAT) 
PAUNGA  CPOSTTIONs VELOCITY } 
[irmoel tPOSslTTONs VELOCITY ) 
Seep NU VECTOR 
NUMSeER OF [TERATIONSs MAX ACCELERATION, GRAVITATION CONSTANT 
Peeve time= STEPs EARTH GRAVITYs FLIGHT TIME 
POWERs SWITCHING NUMBER 

MoRMAT( 3F20.9) 

NITH=W(7s 15) 

V(Bs 20)HW(7s 19)/1090069 

W(Bs 17)=W(7s 16)*W(Rs 20) 
Pew Sree OF [INPUT WOUANTITIES 10 FINAL STORAGE 

W(Bs 16)=2.9%W(Bs 17)/Wl6s 17) : 
SOvMeUm ATION OF EXHAUST VELOCITY 

Pos L=le NIT | 
peeraroON, LIMIT (NIT } 

Mbt CORREC(Ws TT) 

CARE LAMTAB(wHs 7) 

Breer 25 (Wi5s5 I[)s Tal5s 20) 

meeMAlT(3E2028) 

Ge 10 3 

eND 
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Zi 


wn 


10 


I PW 


LAMTAS 

ESGeEL 

loa 

SUBROUTINE LAMTAB(W, T) 


SUBROUTINE COMPUTES, PRINTS bDACKWARD INTEGRATED LAnGDA KiATRIX 


DIMENSION Wl 12s “20s ily ss 8s “GOls Ol Use) 

eee Se rare) 

Ee DS Si 22 aN Sek 

DiGwe Ane lek =) sole 

Rat Otel =. WAC 7 a ie fe ely) 

IGE ae Coa) Se) — Ee ths ome 
PAGS SPACING GF LAMBOA PRINT GUT 

Pah hes 

iS a eae re i, 

Gaelic ane) = Iicoaa, 

aes ee) gt) 

DiOe VO rK =e, 

Otls JI=OC1ls Jy4T( ls Ktls LLIFTC Kae sels Eel) 
INVERTED FINAL VALVE OF FORWARD INTEGRATED LAMBDA iA TR 
MULTIPLIER INTO UAW BDA MATRIX “A 2AGHeshe 

FORMAT( 1°4X« Fl6<6) 

POsmaAl (lows 2a the sO «1 Sale Ae eae ine) 

POG Ad (ClHiUst 1. cs.) (PbO sic 

Or COs © 1 VG.) 

PRINT 2 

PRIN iro se) tls “ler EL ey (ol ls 1) 5 la) 

Pehle, SOs: “ilesc = snr fence = case) 

PRINT 45 A 

ReTURN 

END 
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rc 
3 
xt 


2) 


3 


2 
1 


COPREC 
ASE 
mist 
BSUSROUTINE CORREC(W,s T) 
SUSoROUTINGE COMPUTES NEW INITIAL VALUES FOR NEXT ITERATION 
PIDSeRINDS FINAL VALUES OF INTEGRATED QUANTITIES 
PetENSTON WI12s 20)5 Q(105 10)s Ti?s 8s 50} 
MAE INTEG(Ns T) 
meelO T=ls 6 
WR, 1>)=2e9 
W(45 [+14)=W0 1,5 T)-W(35 1+14) 
Sora ton OF MISS VECTOR AT TARGET (TARGET XP) 
pate>s 1+14)=0.0 
BOyl> J=1. 6 
Pees rhayew( 2, I+1l4)4+wll+5, J)*(Wlls J)—W(35 J+14)} 
Penner ce taRGe! XI] TO EQUIVALENT LAUNCH ERROR (INIT EAE a 
Otis JVENC 1455 J+7) 
Boime TT NVERT (Gs Q) 
eves oLON OF HM STAR MAIRIX 
BO 2°? l=l1s5 § 
pone. J=)]s 6 


Pees It Le jJaew( 5, [F14)—-OQl 1s J)#W(25 J4+14) 


COMPUTATION OF CORRECTED NU VECTOR 
mee tAT(IHISXIBSHINITIAL STATESXIIHFINAL STATE6X12HTARGET STATESX 


Mea - I ALIXOHTARGET XIS8X1LOHINITIAL XI7X1OHNU X 10000) 


pee “47 ( 1A 3 7F 172°) 
Bee er O/ LT ROSXITHISSXI2HMSTAR wiATRIX/IHOELT+e /s6F17-¢6/( lGX6FL76G)) 
peed 5s 3) /26° 
PRINT 3 
ORTHT 2 
Poer> [=15. 29 
XNUFW(55 1)%*10000-0 
PRINT 259 WllsIl)oWl1l, I-14) 9¥(351) 96(351-14) oW0451) oWl25I) »XNU 
PRINT ls AJs (C(wlls J)s1=65 11)s J=8s 13) 
Poeoo |=l1ls 7 
DO 30 J=15 7 
Ras ye CI+5 + 1) 
SEM INVERT Crs Q) 
ear Gege 9941.01 
oO tO T=]. q 
MP9 40 J=2s5 8B 
Miers see =Q0C1ls J-1) 
ys oe y= Te O/W( ds 7} 
TNVERSION OF FINAL VALUE Or FORWARD INTEGRATEU LAMZDA MATRIX 
RETURN 


aus) 


e 
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oO 
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w 


i) 
u) 


kad 


ra 
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\Jl 


ig 
ae 


DPAESG 
aes Se 
eh 
SUBROUTINE INTEG(Ws 7) 
SUSROUTINE CONTKOLS INI TIAE  VAbo SS FOR eG were, 
ANG PRINTS OUT TRAJECTORY VAKT4SbeS 4 baer ee 
PleeRSlOn Wol2s 20s "Dl 1l2s 20) sl eo see 
OY JO J]=1+ 168 
ON As Vi= ee 
ts 1) ae. 
De 20 (| = lease 
Poo sel) = lee 
wes (jee 4 lee + ee) 
G2 se ea og eG) 
Pls 7)=]e6 
e228 c= ene 
Bean, ) =e 
ales 2 Je=sleeO 
Dttis Oa ae) 
BY so. 1 3 rade 
PI ies et OA ee 
PA Ae PeON DLT PONS BOR NT eGR Alero 
Ceol et le yal sea) 
PR 5 VS = a Gier 3s 3s 
jel Sor THRUST EI MTPATIGN (Sw ter Poi 
TGA] 20): eae ee) 
CALL DERV('™. D) 
Poe = (yes) 10 ) 
IMT TAY PGN Or DUNTEGRATT Om Si be esis] 
FORMATCIALIIX4AHTIMec&X1L2nACCELecRATIONTX9OHACCEL MAGBKABHPOS1] TILONDGX 
B4VELOCITYLIX4HMASSITSXSHGAMMA/ZOX9OHX LOCOG/GBEXIHX 10000/G) 
poe Cl ee eae G s “Ors hear 
leer th LA. se Kew ol gore en ee a ere) 
Gg ag pr 
P= ee a 0 eG oO) 
A2=W(os 193/485 20) 
ASHEN(6 9) COV PCBs: ZO) 
ADHD CS (2/7 8658 20 } 
GAME'( 7s 16)/wW(6s 17) 


ev 


Cel 7) | 

Pee =. |e 

Bee Ke Lee 7 

Cs er so ea ees ne) 

Tee 0 ey i ey Opes re © 
PAGE > -fOCING PORT ITRAVECTORY sess 

Poot a 

SO PNT 2 


PRINT 25 W551) sALSAMsW( 151) svills4) oswl1ls7) »GAM 
Pati la Ae lO is 2 ys We de SY eee es ce yt Ls) 
Pe Ce awit Se =H 7s 205) GOs ons 30 
ie oleae oe? ce Nbr FElGHy). Vlse 
Pete (7 6 An oe 1) 
Petey Wy she =e CUO 1) 60s 66) 70 
loo eon BINA STEP LESS “lAANe.O000 I 
CALLM STEP CUWs Ds DEL) 
RETURN 
ell 
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STFP 
LAGEL 
Leo 
SUBROUTINE STEP(wls Dl» DEL) 
SUBROUTINE SELECTs STEP SIZE ON BASIS OF INTEURATION 
ERROR AND THRUST LIMITATION SWITCH POINT 
BPE NSTON J1(240)> D1(240)>5 W2(240)5 D2(240)+ w3'240)> D3(240) 
FLAG=W1(175) 
DEFINITION OF FLAG 
as 2G 
Q1¢1) 
Duleted 
= it 1) 
yaW1(1) 
DEFINITION OF VARIABLES FOR TRIAL INTEGRATION 
Merl PUNKUT(W2,5 D2. DEL) 
TRIAL INTEGRATION 
IF (FLAG—112(186)) 205 155 15 
TEST FOR THRUST LIMIT SWITCH GFF POINT DURING INTEGRATION 
BeAG=0.0 
ee 10 36 
TFC FLAG—“INLF(W20186)—-160¢05 65)) 255 405 40 
meow FOR THRUST LIMIT SWITCH ON POINT DURING INTEGRATION 
Fi AG=1.0 
DFLENDFL¥ASBSF((wW1(187)-Wl 
ee Cllen or STEP < 
wo 35 Jat. 2749 
me Clyso1 tl) 
Peer ya 1 (I) 
Mee 2UNKUT(A2,s D2s DEL) 
INTEGRATION OVER FULL STEP 
Meel RUNKUT(W3> 33. DEL/2.9) 
Meek RUNKUT(W3s 03s DEL/2.0) 
TVO INTFGRATIONS OF HALE STEP FACH 
ee = 92 11 74)7431174)=-160)/ «0009045+.0001 
B(T=S1—2.0) 555 55,5 50 | 
DIFFEPcNCE TEST OF FULL STEP AND HALF STEP INTEGRATION 
DEL=DEL/ST=EST#%.25 
IS (TEST-.95) 605 60, 70 
Pee MIN CwWl(211)> DELJTEST#+#.25) 
SELECTION OF STEP. SIZE FROM ERROR TEST 
mor 75> [=)s 240 
Pe 1 pel) 
Pai t= 37 4 
STORAGE OF INTEGRATED VALUES FOX NEXT INTEGRATION STEP 
Port223)=1TEST#e0090003 
IF (Wi(i75)-FLAG) 80s 90> 80 
W1(212)=W1(175)-FLAG 
W1(175)=FLAG 
STORAGE OF THRUST LIMIT SWITCH POINT 
CALL DERV(W1, D1) 
Pale? 7 =e 6 
RETURN 
END 
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CREEP UERY Cis Dies 
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Wie tee Os 


CAGE Deny (el. 
RETURN 
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Oo CS ee ee aia 


CV Nae 
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el 


eZ 
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20 


25 


BD 


BNVERT 

IPABE L 

lo T 

SUBROUTINE INVERT(Ns QQ) 


SUBROUTINE INVERTS MATRIX BY SIMULTANEOUS DOUBLE PRECISION 
peveREOUCTION OF THE MATRIX TO IDEM AND IDEM TO THE INVERSE 


eee NSION OO(10s5 10), Q@(10, 20) 

Be 10 [=ls 10. 

BG 5 J=l>, 10 

O(lTs J)=QQ(1>s J) 

Q( I» J+10)=0.0 

Qt I» J+#20)=0.0 

Prels J+30)=0.0 

Qt1s 1+10)=1.0 

DO 30 JT=l> N 

DO 14 J=I>5 N 

Mee AOSF(Q( 1s 1))-AS8SF(Q(J>5 1))) lls 14s 14 
Poster oRk LARGEST ELEMENT IN COLUMN 

DO 12 K=ls N 


a= @{ J, ce) 
Qld» K)=QCI]>5 K) 
eels K)=S 


S=Q(J, K+10) 
rs K+10)=O( 1s K+10) 
Q(Is K+10)=S 
TRANSFER ROW OF LARGEST ELEMENT TO FIRST ROW 
CONTINUE 
By=-Q( ls 1) 
DO 15 J=ls N 
Spels J)=OCl, J)/DIV 
mens JtlvUj)=O(ls J+10)/DIV 
DIVISION BY DIAGONAL ELEMENTS 
pO, 3s0 J=]>5 N 
meee) 20s 305 20 
Bv=O(J> 1) 
DO 25 K=ls N 
Os K)=O0(Is K)—$Q01, K)*DIV 
eee lh O)- Ones KEIO )—-OC ls Kt+10)*DIV 
DIAGONALIZATION OF MATRIX 
CONTINUE 
DO 35 T=15-N 
Bo°35 J=ls N 
Seis, J=0Cls J+10) 
RETURN 
END 
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RL=(W(20 L)*#W(25 4)4tW( 259 2)*¥wWl2s Z)4¢Wl2s B)#WC2s 6))#AZ*W( 1» 
GOMANTITY USED IN M STAR MATRIX 

PO 30 I=15 6 . 

DO 39 J=ls5 6 

ey 6 

Doecs K=ls 3 

FHF ( [+59 K+3)%#W(U+559 K+3) 

Pte os J 7)=(F=0(35 1)*#D1(35 J)/D(55 3))#D15s 2)/7A2! 
DERIVATIVE OF M STAR MATRIX 

Pelt S5 J+7)Sswllt+5s J+7)+W( 45 T)FWI4s JV/SRL*¥W(6s 1L7I*WIBs 18) 
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MeOtSs L+7)=W(l+5s9 J+7) 

Peet o> [FP )=Di1+5s J+7) 
Deen wilor So MMe FR] C 

DO 40 T=l>s 168 

Pel y=D(1> 1) 

RETURN 

END 
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DO 2m J=1> 3 
D(1+5s J) ERMZ*#WOT4t5s Jt3)-RL#W(ls J) 
DCT+5s J+3)=—" (1455 J) | 

ADJOINT EQUATIONS (i1ASS INOEPENDENT) 
DBO SO J=)> 6 
worsO J=I1>5 6 
P= Je 
MO. 25 K=l> 3 
FaF+tw(l+59 K+3)*¥W(Ut5s K+3) 
D(I+5, J+7)=(F-D(3» 1)*#0(35 J)*STAR)/GAM 

NERIVATIVE OF M STAR MATRIX 
Besos 1[+4+7)=D(1+5s J+7} 

M STAR IS SYMMETRIC 
Senco [=1s 168 
D1l(I)=DtI». 1) 
PETURN 
END 
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APPENDIX H 
NUMERICAL RESULTS 


H.1 Introduction 


In this appendix the computer output data from a sample run of 
each control mode are reproduced. The data sheets are in pene am 
self explanatory except perhaps for the format of the A. matrix and 


the units. 


All parameters and variables derived in the linear theory are 
printed out. In addition, the important variables are plotted in Pigumes 
Bley oseroyyolel gl oie 


H.2 Format and Units of the Data 
The first pages of each example contain: 
ij tnetime in days irom initiakpom:, 
2), Phe acceleration program in units ot 10°*g., 
3.) The acceleration magnitude in units of 10°“g., 
27) Ue postion ine ua. 
oS.) The velocity in A.u. per day, 
6.) The normalized mass, 
ia) bhevewiltching parameter, ye 
The coordinates are computational coordinates in all cases. 


Following the state and acceleration data are the terminal values 
of interest printed as column vectors. The units are A.u. and 2m 


per day. 


The elements of the M* matrix, which are the next set of numbers, 
are not particularly interesting and may be disregarded. 
Adjacent to the M* matrix is the cost quantity J with units of 


(A.u.)* per (day). 
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The next set of pages contain the elements of the A. matrix for 
both optimal acceleration and optimal thrust programs. For each time 
Point, the@array consists of 50 elements which are interpreted as par- 
tial derivatives or influence coefficients. The array is ordered as in 


the equation 
OS ¢ = SJ (t) cr 


The first 49 elements are the adjoini set for an optimal thrust program. 
mnevadjoint set for an optimal acceleration program is obtained by re- 
placing the seventh column with zeros for the first six elements, and 


meplacing the 1 with the element below it. Thatis: the array 


Mog 9 iste 
Noy A-2og Noe 
T T 


ai 
N. 32 
yields A. (optimal thrust) 


eS ee sees 
No, Neg Ap 
O O 1 
and /\. (optimal acceleration) 
Pi ioe 


ee ope 


T T 
O Oe 


The applicable submatrices, interpreted as partial derivatives, are 
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a 
or, ov, om, 
OV OV OV - 
om om om 


The units are A.u. and A.u. per day for position and velocity, Vie 
seventh column is with respect to a 100% change in mass. For example 


the units of au are A 
om, 100% change in mass 


The first set of data is for unrestricted VSI control. The second 


set is for CSI control with a = tex lee 
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Fig. H-a. Acceleration schedule for 150-day Earth-Mars transfer. 
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Fig. H-d. Sensitivity of r, (t¢) to position variations. 
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Fig. H-e. Sensitivity of r, (t¢) to velocity variations. 
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Fig. H-f. Sensitivity of ry (t¢) to position voriations. 


io2 





140 


ISO 


hoe Ory (te) (METERS) 

24 Ov, (METERS/SECOND) 
rf Ory (t¢) 

25 O vy 
Ory (t ¢) 


Ovz 


hog = 





0 10 20 30 40 50 60 70 80 90 100 !10 120 130 140 150 
TIME (DAYS) 


Fig. H-g. Sensitivity of ry (tg) to velocity variations. 
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Fig. H-h. Sensitivity of r, (t¢) to position variations. 
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Fig. H-i. Sensitivity of r, (ts) to velocity variations. 
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Fig. H-j. Sensitivity of v, (t¢) to position variations. 
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Fig. H-k. Sensitivity of v, (t¢) to velocity variations. 
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Fig. H-l. Sensitivity of Vy (t¢) to position variations. 
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Fig. H-m. Sensitivity of Yq (t¢) to velocity variations. 
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Fig. H-n. Sensitivity of v, (ts) to position variations. 
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Fig. H-o. Sensitivity of v, (tf) to velocity variations. 
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Fig. H-p. Sensitivity of r (t¢) to mass variation for optimal thrust program. 
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